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DYNAMIC  STRUCTURAL  DATA  ANALYSIS 


1 .  INTRODUCTION 

This  document  describes  methods  and  procedures  that  deal  with  the  ^  ■ 
reduction,  presentation,  and  analysis  of  environmental  data  which  have  speci¬ 
fic  application  to  structural  evaluations  and  fall  under  the  following 
categories:  ■ 

a.  Vibration 

b .  Shock 

c.  Acoustics 

d.  Strain 

Engineers  and  other  personnel  involved  in  structural  data  analysis 
are  aware  that  no  one  method  of  reduction  or  presentation  is  always  the  best 
one.  Usually,  the  correct  method  depends  on  the  type  of  data  and  the  use  to 
be  made  of  the  results. 

The  data  analyst,  generally,  has  the  necessary  equipment  at  his 
disposal  to  obtain  the  data  he  desires.  However,  the  equipment  limitations 
and  characteristics  must  be  weighed  carefully  against  conflicting  (usually 
mathematical)  requirements.  For  this  reason,  the  mathematics  involved  in  the 
analysis  are  integrated  into  the  operational  discussion  as  much  as  possible. 

In  any  type  of  data  analysis,  knowledge  of  the  mathematical  formulae  which 
govern  the  processes  involved  will  produce  a  clearer  insight  into  the  analysis 
(or  analyses)  being  conducted. 

Basically,  this  document  deals  only  with  the  essential  background 
material  and  methods  used  in  the  analysis  of  a  prerecorded  signal.  It  is 
assumed  that  the  results  of  the  analysis,  whether  it  be  spectral  or  statistical, 
will  be  passed  on  to  a  cognizant  design  or  test  engineer. 

2.  VIBRATION  DATA  ANALYSIS 

Vibration  data  may  be  measured,  reduced,  and  presented  in  two 
different  ways.  One  method  is  to  consider  the  magnitude  of  the  vibration  at 
a  given  frequency  of  vibration.  In  other  words  construct  a  plot  of  magnitude 
versus  frequency.  The  other  type  of  presentation  is  one  of  magnitude  of 
vibration  at  a  given  frequency  over  a  period  of  time,  i.e.  magnitude  versus 
time  plots  as  shown  in  Figure  1. 

A  plot  of  both  time  and  frequency  versus  magnitude  is  shown  in 
Figure  2,  and  illustrates  the  difference  between  the  two. 
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Figure  1.  Typical  Vibration  Time-History  Broadband  Random  Vibration. 
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Figure  2.  Vibration  Magnitude  Surface. 
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2.1  TYPES  OF  FREQUENCY -MAGNITUDE  PLOTS 

A  presentation  of  the  magnitude  of  a  vibration  versus  the  frequency 
of  a  vibration  is  called  a  frequency  spectrum.  The  analysis  of  a  periodic 
wave,  as  shown  in  A  of  Figure  3  when  performed  over  a  large  finite  number  of 
cycles,  will  produce  a  line  spectrum  or  discrete  frequency  spectrum  as  shown 
in  B  of  Figure  3.  Instead  of  a  graphical  presentation,  the  results  of  this 
type  of  analysis  might  well  be  presented  as  a  tabulation  of  frequencies  and 
associated  amplitudes . 


A.  Typical  Vibration  Time-History  Sum  of  Two  Sinusoids 


FREQUENCY 


B.  Results  of  Data  Analysis  of  Time-History  of  A.,  Sum  of  Two  Sinusoids 


Figure  3.  Typical  Vibration  Time-History  and  Data  Analysis. 


This  type  of  presentation  might  represent  vibration  magnitudes 
through  fixed  filter  bandwidths  at  an  instant  in  time.  For  example,  this 
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could  be  accomplished  by  passing  the  instantaneous  acceleration  signal 
through  a  set  of  parallel  filters  and  recording  (on  an  oscillograph)  each 
filter  output  simultaneously.  The  data  required  for  the  frequency  spectrum 
at  time  ti,  would  then  be  taken  by  observing  on  each  oscillograph  trace  the 
amplitude  present  at  t]^.  (See  Figure  4.) 


Sometimes  it  is  desired  that  the  amplitude  associated  with  the  fre¬ 
quency  spectrum  be  a  time-averaged  value.  For  example,  it  might  be  desirable 
to  plot  the  root  mean  square  (rms)  acceleration  experienced  at  a  particular 
station  of  a  missile  during  the  free  flight  portion  of  a  given  missile  maneuver. 
In  this  case,  it  would  be  necessary  to  average  the  instantaneous  acceleration 
in  some  way,  over  the  time  of  interest.  In  some  cases,  obtaining  data  related 
to  an  instant  in  time  can  be  precluded  by  averaging  time  constants  and  damping 
built  into  the  instrumentation  being  used. 

Although  the  presentation  form  discussed  previously,  which  involves 
frequency  as  the  abcissa,  can  be  referred  to  as  "spectral  analysis",  the 
concept  most  often  associated  with  this  term  is  that  of  the  "power  spectrum". 

In  general,  vibration  data  are  aperiodic  (i.e.,  random  in  nature). 
Such  a  vibration  signal  cannot  be  resolved  into  discrete  sinusoidal  com¬ 
ponents  such  as  amplitude  and  phase  spectra.  However,  the  energy  associated 
with  a  random  time  function  taken  over  a  finite  length  of  time  is  finite. 

Thus,  many  data  can  be  obtained  for  an  arbitrary  signal,  if  the  distribution 
of  its  power  with  respect  to  frequency  is  known. 

Actually,  it  is  the  rate  of  change  of  average  power  with  frequency 
which  is  of  prime  interest  and  this  introduces  a  presentation  concept  of 
"power  spectral  density". 
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2.2  MATHEMATICAL  RELATIONSHIPS  THAT  DESCRIBE  VIBRATION  PHENOMENA 

2.2.1  Sinusoidal  Functions 

The  simplest  time-varying  quantity  of  concern  in  vibration  is  a 
sinusoidal  function  of  constant  peak  amplitude.  Am  and  frequency,  f. 

a(t)  =  Am  sin  (tut  +  0)  (1) 


where  u)  =  2TTf 

0  =  phase  angle  with  respect  to  the  time  origin. 

2.2.2  Periodic  Functions 

A  function  which  has  a  waveform  which  repeats  itself  exactly  during 
each  time  period  t  such  that : 


a(t)  =  a(t  +  nr)  (2) 

where  n  is  any  integer,  then  it  is  said  to  be  a  periodic  function  having  a 
fundamental  frequency  where  f  equals  l/j.  A  periodic  function  may  be  repre¬ 
sented  by  the  line  spectrum  as  shown  in  Figure  5. 


Figure  5.  Line  Spectrum  to  Represent  Fourier  Series  of  Periodic  Function. 

Fundamental  Frequency  fg  =  1/t,  where  r  is  Period  of  Function. 

2.2.3  Phase  -  Coherent  Functions 

A  phase -coherent  function  is  defined  as  a  sine  wave  whose  amplitude 
and  frequency  vary  smoothly  with  time.  It  is  assumed  that  the  changes  in 
amplitude  and  frequency  are  significant  in  magnitude  only  over  a  time  duration 
equal  to  a  number  of  cycles  of  the  waveform.  It  is  expressed  as: 

a(t)  =  ao(t)  sin  [tiu(t)  +  0]  (3) 

where  and  w  are  functions  of  time  and  9  is  the  phase  angle  of  the  waveform 
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at  the  origin. 

In  a  strict  sense,  phase -coherence  is  determined  by  a  single 
reference  frequency  which  persists  for  all  time.  This  type  of  function  is 
useful  in  describing  a  vibration  measured  in  a  system  whose  material  fre¬ 
quency  varies  with  time. 

2.2.4  Complex  Functions 


A  complex  function  is  a  time-history  which  consists  of  a  sxnn  of 
sinusoidal  functions  but  which  is  not  necessarily  periodic.  It  can  be  repre¬ 
sented  graphically  by  a  line  spectrum  such  as  that  shown  in  Figure  5.  The 
lines  are  not  equally  spaced  along  the  frequency  axis  for  a  complex  function, 
and  if  the  function  consists  of  a  continuous  distribution  of  sinusoidal 
functions,  it  is  best  described  by  the  Fourier  spectrum. 

2.2.5  Random  Functions 


A  random  time  function  is  composed  of  a  distribution  of  sine  waves 
which  is  continuous  and  contains  all  frequencies.  It  cannot  be  represented 
in  a  deterministic  sense;  that  is,  if  its  value  at  a  given  instant  of  time  is 
known,  it  is  known  only  that  at  some  time  A  t  seconds  later  there  is  a  certain 
probability  that  the  value  will  be  within  a  specific  range  of  values.  Random 
functions  can  only  be  described  in  statistical  terms  as  opposed  to  a  sine  wave 
which  is  completely  deterministic. 

2.2.6  Transient  Function 

A  transient  function  is  a  time -varying  function  which  is  nonzero 
only  over  a  given  restricted  time  interval.  Shock  is  a  force  or  acceleration 
represented  by  a  discontinuous  function  of  time.  Physically,  no  force  or 
acceleration  is  strictly  discontinuous  and,  therefore,  shock  also  refers  to 
those  forces  and  accelerations  that  have  first  time  derivatives  of  very  high 
magnitude  and  are  not  sinusoidal. 

2.2.7  Fundamental  Statistical  Formulas 


a.  Mean  Value.  The  mean  value  of  a  function  x(t)  is  defined  by: 

T 

Ic  =-^J  x(t)dt  (4) 

o 

Where  T  is  the  time  interval  over  which  the  data  points  x(t)  are  taken.  In 
most  vibration  cases  the  mean  value  will  be  zero. 

b.  Average  value.  The  term  "average  value"  of  a  function  x(t) 
is  employed  to  denote  the  time  average  of  the  magnitude  of  a  function  with 
zero  mean  value: 
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c.  Mean  square  value, 
is  the  time  average  of  the  square 


The  mean  square  value  of  a  function  x(t) 
of  the  function: 


(6) 


d .  Root  mean  square  value  (RMS).  The  RMS  value  is  the  square 
root  of  the  mean  square  value. 

2.3  Analysis  of  Periodic  Functions 

2.3.1  Fourier  Series 

Any  periodic  function  f(0),  with  period  2rr,  which  satisfies  the 
Dirichlet  conditions  may  be  expanded  into  a  Fourier  series 


<^0  r 

f(0)  =  'Y~  ^n 

n=l 


(7) 


which  consists  of  a  sum  of  sine  waves  whose  frequencies  are  all  multiples  of 
the  fundamental  frequencies.  The  Dirichlet  conditions  to  be  met  by  the 
function  are: 

a.  It  must  have,  at  most,  a  finite  number  of  discontinuities  in 
one  period. 

b.  It  must  have,  at  most,  a  finite  number  of  maxima  and  minima  in 
one  period. 

c.  The  integral  J  |f(0)j  dS  finite. 

-TT 

The  coefficients  a^j  and  b^j  are  given  by: 


^n 


f(9)  cos  n9  d0 


(8) 


bn 


sin  n9  d9 


(9) 


The  series  can  be  expressed 


alternatively  as: 
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% 

00 

ao 

f(t)  = 

n=l 


where  C^j 
il-n 


^  +  bji^  =  the  Fourier  coefficients 

^n  -x 

tan  ( -r —  1=  the  relative  phase 
'  ^n  y 


(11) 

(12) 


In  each  case,  Sq  is  equal  to  the  average,  or  direct  current  (d-c),  value  of 

T~ 

the  function  under  study,  A  plot  of  jCj^j  versus  ntr  such  as  shown  in  Figure  6 
is  called  the  frequency  spectrum  of  the  given  f(t).  For  a  periodic  function, 
this  graph  is  a  discrete  frequency  spectrum,  as  opposed  to  a  continuous 
spectrum  such  as  arises  in  the  study  of  aperiodic  signals, 

2 


Figure  6.  Frequency  Spectrum  of  f(t)  in  Figure  7. 


i/ 


If  the  given  periodic  function  exhibits  certain  symmetry  conditions, 
several  simplifications  can  be  applied  to  the  analysis. 

If  the  function  is  even,  such  that 

f(t)  =  f(-t),  (13) 

then  the  Fourier  series  expansion  contains  only  cosine  terms  (b^^  =  0)  and  the 
constant  term  (if  a^  ^  0) . 

If  the  function  is  odd,  such  that 

f(t)  =  -f(-t),  (14) 


I 
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then  the  Fourier  series  expansion  contains  only  sine  terms  (aj^  =  0) . 

The  periodic  function  f(9)  with  period  2tt  contains  only  even 
harmonics  if  it  satisfies  the  condition  that: 

f(e  ±  tt)  =  f(e)  (15) 

Likewise,  the  periodic  function  contains  only  odd  harmonics,  if  it  satisfies 
the  condition  that 


f(6  ±  tt)  =  -f(e)  (16) 

For  example:  Find  the  Fourier  series  expansion  of  the  periodic  rectangular 
wave  of  Figure  7: 


Figure  7.  A  Periodic  Rectangular  Wave. 


which  is  expressed  over  one  period  ^  <  t  ^  • 


T  T 

0  ~<  t  <  — 


f(t)  =  —  <t<  ^ 

„  T  T 

0  <  t  <  Y 


Using  equation  (8),  we  determine  the  a^,  coefficient  for  cosine  terms  as; 


(17) 
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(18) 


(19) 


(20) 


(21) 


By  inspection  of  the  even  symmetry  of  the  given  wave,  it  can  be 
seen  that  the  =  0  for  this  case. 


The  average  value  can  be  determined  by  substitution  of  n  =  0  in 
=  2A/nTT  sin  (nTT/2),  which  then  takes  the  indeterminate  form  o/o.  This  can  be 
evaluated  using  L'Hospital's  rule  or  the  average  value  in  this  case  can  be 
determined  to  be_^  .  The  total  series  can  then  be  expressed  as: 

2 


f(t)  = —  (cos  cut  -  cos  3uut  +  cos  5u;t 
2  TT  J  5 


Fourier  series  can  also  be  written  in  exponential  form  as: 

"hx> 


...) 


(22) 


£<0  =  I 


jncot 


one 

n=-oo 


(23) 


where 


I. 

.2 


On 


-I 

2 


with  n  ranging  over  all  values  from  -«>  to  +». 
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An  advantage  of  the  exponential  form  of  the  Fourier  series  is  that  is  provides 
a  more  convenient  transition  to  Fourier  integrals  and  transforms. 

2.3.2  Harmonic  Analysis 

Harmonic  Analysis  is  the  determination  of  the  coefficients  of  the 
series  expansion  for  periodic  waves  given  in  the  previous  section.  One 
method  of  harmonic  analysis  uses  parallel  bandpass  filters.  The  recorded 
signal  is  played  back  through  each  filter  in  turn,  and  the  filter  output  is 
again  recorded  for  further  analysis  on  an  instrument  such  as  a  level  recorder 
or  an  oscillograph.  Another  method  widely  used  takes  advantage  of  the 
"heterodyning"  effects  of  a  variable  oscillator  as  it  is  mixed  with  the  input 
data,  and  then  filtered.  The  filtering  is  performed  at  either  a  fixed  inter¬ 
mediate  frequency  or  at  a  zero  IF.  The  local  oscillator  is  made  to  sweep 
automatically  at  some  preselected  rate.  A  tape  loop  input  is  always  used  with 
this  type  of  analyzer,  whereas,  in  the  parallel  filter  concept,  a  simple 
tape  input  scheme  is  acceptable. 

2.3.3  Statistical  Analysis 


Statistical  analysis  of  shock  and  vibration  data  is  utilized  to 
determine  the  "amplitude  spectrum"  of  stochastic  or  random  functions.  It 
largely  concerns  the  predetection  or  probability  that  the  vibration  data  will 
fall  between  arbitrary  levels.  From  this  analysis  a  distribution  function  is 
usually  derived  that  will  permit  an  evaluation  of  vibration  and  shock  effects. 

2. 3. 3.1  Density  Functions 

In  general,  it  is  desirable  to  produce,  in  graphic  form,  a  pre¬ 
sentation  from  which  can  be  obtained  a  discrete  probability  associated  with 
a  certain  event.  (For  example,  the  probability  of  occurrence  of  a  certain 
range  of  accelerations.)  One  concept  which  can  be  introduced  to  do  this 
is  that  of  the  probability  density  function. 

By  way  of  approaching  the  concept  of  probability  density,  assume 
that  the  distribution  of  height  among  American  males  must  be  determined. 

The  height  of  a  man  may  take  on  any  value  within  a  specified  interval  and 
thus  represents  a  continuous  random  variable.  (As  also  does  most  inflight 
vibration  data.)  Assume  that  a  representative  sample  of  1000  men  is  selected 
for  this  determination.  Then,  their  heights  are  measured  and  they  are 
grouped  according  to  the  nearest  even  inch  (5  feet,  0  inches;  5  feet,  2  inches 
5  feet,  4  inches;  etc.).  All  heights  from  4  feet,  11  inches,  to  5  feet,  1 
inch,  will,  for  example,  be  grouped  in  the  5-foot,  0-inch  category.  The 
relative  frequency  of  men  found  in  each  height  interval  is  then  the  number 
grouped  in  that  interval  divided  by  the  total  number  (1000) .  A  typical 
plot  of  such  a  height  distribution  is  shown  in  Figure  8.  In  Figure  8,  x 
represents  the  height,  n^/n  the  relative  number  of  men  grouped  in  an  interval 
A  X  inches  about  x  as  the  mean,  n  is  the  1000-man  sample.  Since  all  heights 
in  a  2 -inch  interval  about  the  even  heights  are  grouped  together,  they  are 
shown  as  horizontal  lines  covering  the  2-inch  grouping.  For  a  large  enough 
sample,  it  may  be  stated  that  n^/n  represents  the  probability  P(xj),  that 
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(HEIGHT) 

Figure  8.  Relative  Distribution  of  Heights. 


the  height  of  an  American  male  will  be  between: 


Ax  ,  .Ax 
_  and  Xi  +  __ 
o  J  2 


To  do  away  with  the  dependence  of  the  size  of  the  interval  Ax  chosen 
for  the  height  distribution,  the  height  distribution  is  represented  at  a 
particular  height  and  of  interval  Ax  by  a  rectangle  of  area  n^/n.  The  height 

of  this  rectangle  will  be  width  will  be  Ax. 

Such  a  relative -frequency  curve  is  called  a  histogram.  Two  histograms  for 
the  height  example  used  here  are  shown  in  Figure  9  (where  Ax  equals  two 
inches)  and  Figure  10  (where  Ax  equals  one  inch). 
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For  a  large  enough  sample  size,  the  probability  that  a  man's  height  will  lie 
between  two  given  heights,  (e.g.,  3  feet  7  inches,  and  5  feet,  9  inches) 
will  now  be  the  area  of  the  histogram  between  these  two  limits. 


Figure  10.  Histogram  of  Height  Distribution,  Ax  Equals  One  Inch. 


If  the  histogram  approaches  a  smooth  curve  as  Ax  approaches  zero, 
the  ordinate  takes  on  the  form  of  a  probability-density  function,  with  the 
area  under  the  curve  between  two  points  giving  the  probability  that  the 
height  will  be  found  between  these  two  points.  The  probability-density 
function  as  the  limiting  case  of  the  two  histograms  of  Figures  9  and  10,  is 
shown  in  Figure  11. 

The  area  between  the  two  points,  5  feet,  6  inches,  and  5  feet,  8 
inches,  shown  crosshatched  in  the  figure,  represents  the  probability  that  a 
man's  height  will  be  found  in  that  range. 


Figure  11.  Probability-Density  Function  Corresponding  to  Height  Histogram. 
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Thus,  what  has  been  established  then,  is  a  function  p(x)  where: 


p(x)  = 


lim 

Ax"K) 

n-te> 


nx/  tx 
n 


(25) 


where  p(x)  is  the  probability-density  function.  Px  represents  the  number  of 
samples  of  the  total  n  found  in  the  range  x  -  A  x/2  to  x  +  Ax/2.  To  then  find 
the  probability  that  a  value  is  in  a  specified  range,  x^  to  xg ,  p(x)  is 
integrated  through  this  range  by: 


Prob  (x^  <  X  <  Xg)  =  J  p(x)  dx 


(26) 


Since  it  is  certain  that  every  measurement  must  yield  some  real  value,  it  is 
necessary  to  have: 

J  p(x)  dx  =  1  (27) 

-OD 


The  histogram  or  probability-density  curves  must  be  normalized  to 
have  unity  area. 

To  determine  the  statistical  averages  from  the  probability-density 
function  just  described,  it  can  be  reasoned  that  in  the  limit,  as  the  number 
of  measured  values  becomes  large  and  as  Ax  approaches  0,  then: 


av 


F(x)  -  J 


F(x)  p(x)  dx 


(28) 


where  the  abbreviation  "av  F(x)"  is  used  to  mean  the  average  value  of  F(x) 
approached  by  a  large  number  of  trials.  Other  designations  in  use  include 
"<  F(x)  >"  and  "mathematical  expectation  of  F(x),  written  E[F(x)]. 

As  a  special  case,  the  average  value  of  x  is: 


00 

av  X  =  J  X  p(x)  dx  (29) 

-00 


The  average  value  of  x,  av  x ,  is  also  called  the  first  moment,  m^ ,  of  x  by 
analogy  with  the  concept  of  moments  in  mechanics.  In  mechanics,  the  first 
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moment  of  a  group  of  masses  is  just  the  average  location  of  the  masses  or 
their  center  of  gravity.  For  example,  the  center  of  gravity  of  the  masses 
m^ ,  ma  ...  mg  ,  in  Figure  12  can  be  found  simply  by  taking  amounts  about  the 
point  where  x  equals  0.  The  average  value  of  x,  or  the  center  of  gravity, 
is : 


av  X 


m^  Xj  +  nts  xa  •  •  • 
m^  d”  mg  ... 


(30) 


If  the  mass  is  not  concentrated  as  discrete  points,  but  is  continuous  over 
the  bar,  the  center  of  gravity  can  be  calculated  by  considering  a  differential 
mass  m(x)  dx  [with  m(x)  the  linear  mass  density]  located  x  units  from  the 
origin.  The  center  of  gravity  is  then  found  by  summing  all  the  mass  contri¬ 
butions.  This  sum  becomes  an  integral  to  yield: 

1 

m^  =  av  X  =  J  xm(x)  dx  (31) 

0 


The  second  moment  in  mechanics  is  just  the  moment  of  inertia  of  a  mass  or  the 
turning  moment  of  a  torque  about  a  specified  point.  By  analogy  with  mechanics 
the  second  moment  of  a  random  variable  is  the  average  value  of  the  square  of 
the  variable.  For  a  continuous  variable,  this  is  given  by: 


CO 

mg  =  av  x*  =  J  X*  p(x)  dx  (32) 

—00 


Higher  moments  can  be  defined,  but  in  general,  the  nth  moment  of  x  or  av  x 
is  given  by: 


CO 

nijj  =  av  x’^  =  x’^  p(x)  dx 

-GO 


(33) 


The  density  function  p(x)  plays  the  role  of  a  weighting  function  throughout. 

One  significance  of  these  moments  in  voltage  measurements  for  in¬ 
stance  is  that  mj  gives  the  d-c  value  whereas,  mg,  is  found  to  give  the  mean 
squared  voltage  or  the  mean  power. 

2. 3. 3. 2  Distribution  Functions 

/ 
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A  second  important  function  called  the  distribution  function  p(x), 
is  defined  as  the  probability  that  the  value  is  less  than  some  specified  x. 

It  is  given  by: 

X 

p(x)  =  [  p(x)  dx  (34) 

“00 


From  the  non-negative  character  of  the  density  function  p(x),  it  can  be  seen 
that  p(x)  cannot  decrease  with  increasing  x  and  also  that 

p(-oo)  =  0  and  p(“)  ~  1 


p(x)  is  a  continuously  or  monotonically  increasing  function  going  from  0  to  1 
as  shown  in  Figure  13.  At  any  point  where  the  random  variable  has  a  finite 
probability  the  function  p(x)  has  a  discontinuity  equal  to  the  amount  of  the 
respective  probability.  The  value  of  the  distribution  function 


/at  the  minimum  possible  value  of  x  is  p(^ JLiHill?)  =  1.0  since  x  is 
z  y  a  ■’ 

never  less  than  x  min;  for  the  maximum  possible  value  of  x,  p  is 


zero  since  x  never  exceeds  x  max.  Since  the  probability  of  exceeding  that 
value  is  one  minus  the  probability  of  not  exceeding  that  value,  the  ordinate 

of  Figure  19  may  be  changed  to  p  *^-^if  the  scale  is  changed  to  vary  from  1 

to  0  instead  of  from  0  to  1. 


M+Mg+M^+M^+Mg  =  I 


Figure  12.  Discrete  Mass  Distribution. 


The  term  a  (sigma)  represents  the  standard  deviation  of  the  random 
variable  from  the  mean  value  of  the  frequency  function.  It  is  used  as  a 
measure  of  the  spread  of  the  density  function  about  the  mean  value,  m,  and  is 
the  square  root  of  the  variance  or  second  central  moment  ye  • 
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2. 3. 3. 3  Comparison  of  Sinusoidal  and  Random  Functions 

Although  a  sinusoidal  function  is  deterministic,  the  probability 
density  and  probability  distribution  of  a  sine  wave  may  be  determined  for 
comparison  with  those  of  a  random  function.  The  conparison  may  be  made  on 
the  basis  of  (1),  the  instantaneous  values  x|  tj  of  the  function,  or  (2)  the 
peak  values  or  maxima  Xp(t)  of  the  function.  It  will  be  seen  in  the  following 
paragraphs  that  although  the  instantaneous  values  of  a  truly  random  wave  can 
approach  infinity,  the  probability  density  function  of  the  random  wave  under 
consideration  is  well  behaved  throughout  the  entire  region  of  observation. 

On  the  other  hand,  although  the  instantaneous  values  of  a  sine  wave  are  well 
behaved  and  in  fact,  deterministic,  the  probability  density  function  becomes 
asymptotic  to  infinity  at  the  peak  amplitudes  of  the  sine  wave.  An  investi¬ 
gation  would  reveal  that  although  the  distribution  function  in  all  cases  is 
well  behaved,  several  probability  density  functions  display  discontinuities. 

In  general,  as  the  crest  factor  (ratio  of  peak  to  rms  value)  of  a  signal 
decreases,  Jthe  probability  density  function  begins  to  increase  indicating 
the  occurrence  of  one  or  more  discontinuities. 


Figure  13.  Probability-Distribution  Function. 


2. 3. 3. 4  Distribution  of  Instantaneous  Values 

A  comparison  of  some  of  the  identifying  characteristics  of  the 
probability  densities  and  distributions  for  the  instantaneous  values  of  a 
sinusoidal  function  and  the  particular  case  of  Gaussian  random  vibration 
is  shown  in  Figures  14  through  17.  The  term  Gaussian  is  used  to  describe 
a  random  function  whose  instantaneous  value  is  defined  by  a  probability 
density  which  is  proportional  to  the  exponential  of  a  negative  quadratic 
function  of  the  values  of  the  variable.  It  is  one  of  the  most  important 
distributions  in  noise  theory  and  is  expressed  for  one  variable  as: 
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\  _  1  -(x-x 

P(x)  =  -  e  o 

a  /Tn 


(35) 


Here  the  parameters  have  been  adjusted  to  give  the  required  condition: 


00 

["  p(x)  dx  =  1 
(• 

-00 


(36) 


and  also: 

00 

=  I*  X  p(x)  dx  =  x^  (37) 

-00 

00 

Ola  =  I*  x^  p(x)  dx  =  X  ^  +  o® 

-00 

00 

lie  =  ^  (x-TOj  )*  p(x)  dx  =  ii^  - 

-oo 

That  is,  x^  is  the  mean,  o*  is  the  variance,  and  a  is  the  standard  deviation. 

An  important  theorem  in  statistics  called  the  central  limit  theorem 
shows  under  general  assumptions  that  the  distribution  of  the  sum  of  an 
indefinitely  large  number  of  other  independently  distributed  quantities  must 
approach  the  Gaussian  distribution  no  matter  what  the  individual  distributions 
may  be.  In  other  words,  if  x^  ,  Xj ,  ...  x^^,  are  independent  random  variables 
with  arbitrary  distributions,  the  sum: 


X^  +  Xg  +  ...  Xjj 

X  =  -  (40) 

vAT 

approaches  a  normally  distributed  variable  as  n  approaches  infinity.  If  Xj^ 
has  mean  zero  and  variance  of  <  <»,  then  x  has  mean  zero  and  variance 
0^  where : 


0^  = 


aj®  +  +  ...  +  On* 


(41) 
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# 

The  probability  density  function  of  the  instantaneous  values  of  a 
random  variable  following  the  normal  or  Gaussian  distribution  is  shown  by 
curve  A  in  Figure  18.  To  express  the  probability  density  function  of  the 
instantaneous  value  of  a  sinusoid,  several  ideas  in  functional  relationships 
must  be  introduced. 

In  general,  if  the  probability  density  function  q(6),  of  a  variable 
6  is  available  and  it  is  desirable  to  calculate  the  probability  density 
'  function  p(x)  for  x  equals  F(6),  it  is  convenient  to  invert  the  functional 

relationship  and  write  0  =  f(x).  Then  by  direct  substitution: 

q(e)d6  =  q[f(x)]f'(x)  dx  (42) 

where  f  (x)  = - ; - . 

dx 


-A  A 

Figure  14.  Probability  Density  Function  for  Sinusoidal  Distribution. 
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Hence  if  f  (x)  is  a  single-valued  function  of  x,  then: 

p(x)  =  q[f (x)] If '(x) 1  (43) 

If  f(x)  is  multiple  valued,  the  expression  on  the  right  is  summed  over  all 
the  values.  For  the  sinusoidal  case  let 

X  =  A  sin  6  (44) 

where  0  has  a  uniform  distribution  in  range  to  .  Then  the  values  of 

X  are  distributed  like  those  of  the  ordinates  of  a  sine  wave  with  amplitude  A. 
To  find  the  probability  density  function  for  x: 


q(0)  = 


0  when  0  <  and  0  >  ^ 


1  ,  ■■n  ^  ^  3n 

- —  when  —Tr<Q<  — 
2tt  2  2 


where  ensures  a  normalized  condition . 
2tt 


Then  solving  equation  (44)  for  q  as  a  function  of  x: 


0  =  f(x)  =  arc  sin  — 

A 


f'(x)  =  =  [A^-X^*]"^ 


Note  that  each  value  of  x  in  the  range  -A  to  A  corresponds  to  two  values  of  0 
and  hence  that  the  right  side  of  equation  (43)  must  be  multiplied  by  two  to 
obtain  p(x).  Then: 


p(x)  = 


0  when  x  <  -A  and  x  >  A 


-..2  \  -h 


-x= ) 


when  -A  <  X  <  A 
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A  plot  of  p(x)  =  "  is  shown  by  curve  B  in  Figure  18  and  also  in 

Va“  -  X®  tt 

Figure  14.  Note  that  at  x  =  +  A,  p(x)  blows  up.  The  average  value  is  given 
by: 

^  xdx 

-  J  wp-r?  ■  “  («) 

-A 


and  the  mean  square  is : 


tns 


The  mean  square  value  could  also  have  been  obtained  from  the  original 
function  given  that  it  was  periodic  with  period  2n. 

T 

|xl^  "  Y  J  I’  (51) 

“  0 
2rr 

1x1“  J  1a  sin  el®  de  =  ^  (52) 

0 

Likewise  the  distribution  function  is: 


fT  ^ 


arc  sin- 


Thus,  it  can  be  seen  that  it  is  possible  for  the  probability  density  to  be¬ 
come  infinite  provided  the  integral  over  a  finite  interval,  including  the 
singularity,  exists. 

The  probability  density  functions  of  equations  (35)  and  (48)  are 
symmetrical  about  their  means.  If  the  mean  is  assumed  to  be  zero,  then  the 
probability  that  x  exceeds  a  given  absolute  value  (or  magnitude)  1x1  is  twice 
the  probability  that  it  exceeds  the  same  absolute  value  in  either  the  positive 
or  negative  sense.  Therefore,  it  is  convenient  to  plot  the  probability 
distribution  function  in  terms  of  the  absolute  value  of  x,  i.e.. 
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p  (|xl  >)  or  p  1^  >  ^ 


as  shown  in  Figure  19 . 


The  notation  —  or  —  is  often  used  for  the  abscissa  value  in  place 
a  a 

of  X  or  Ixj.  Division  by  the  standard  deviation  in  a  sense  "normalizes"  the 
plot  in  as  much  as  it  reduces  the  abscissa  value  to  a  dimensionless  quantity 
which  holds  its  significance  no  matter  what  the  value  of  the  random  variable 
The  probability  distribution  functions  for  the  Gaussian  and  sinusoidal 
functions  are  obtained  by  integration  of  equations  (35)  and  (48),  and  are: 


Gaussian 


-(xiV2o») 

€ 


a 


(54) 


where  x^  is  a  dummy  variable. 


Sinusoid 


< 


a 


-  y 


arc  cos 


(y/F 


(55) 


The  probability  associated  with  the  Gaussian  distribution  function  above  is 
essentially  twice  the  area  under  either  side  of  the  normal  curve  from  the 
point  in  question  out  to  the  extreme  value  of  the  respective  side.  It  could 
also  be  expressed  as 


x/cr 

(1  ;)■ 

0 

(56) 


The  probability  that  —  is  greater  than  a  given  normalized  value  is  written 

conversely  the  probability  that  —  is  less  than  a  given  normalized 


value  is  written  Also,  p(^-^  ^  ^  =  1  -  •  The  probability 

that  an  instantaneous  voltage  which  follows  the  normal  distribution  will  be 
below  some  value  Ka  (where  K  is  a  constant)  can  be  expressed  in  terms  of  the 
error  function  as: 
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K 


Prob  ( jx 1  <  Ka)  =  ^ 


(57) 


The  error  function  is  described  mathematically  as; 

2  ^ 


erf  X  =  — 

Vn 


-y® 

e  ^  dy 


(58) 


and  is  tabulated  in  many  books  on  statistics  and  in  mathematical  tables. 

The  cumulative  probability  distribution  can  be  related  to  the  error 
function  as  follows; 


P(x)  =  h 


1  +  erf 


(x-a) 

ca/F 


(59) 


where  p(x)  is  the  cumulative  distribution  and  a  is  the  mean  of  the  normal 
distribution  under  study. 

2. 3. 3. 5  Distribution  of  Peak  Values 

When  the  peak  values  or  maxima  of  a  function  are  considered,  the 
two  statistical  functions  differ  from  those  found  for  the  instantaneous 
values.  For  a  sine  wave,  all  maxima  are  of  equal  magnitude  and  the  probability 


density  function  p(Xp/a)  becomes  a  Dirac  delta  function  as  shown  by  curve  B 


in  Figure  20.  For  broadband  Gaussian  vibration,  i.e.,  vibration  with  non¬ 
zero  spectral  density  over  a  frequency  bandwidth  which  is  not  small  compared 
to  the  average  or  center  frequency  of  the  bandwidth,  the  distribution  of 
peak  values  is  normal  as  shown  by  curves  A  in  Figure  18  and  by  the  curve  in 
Figure  16.  However,  for  narrow  band  Gaussian  noise,  i.e.,  noise  with 
negligible  spectral  density  except  in  a  frequency  bandwidth  which  is  small 
compared  to  the  center  frequency,  the  distribution  of  peak  values  becomes  the 
Rayleigh  distribution.  The  probability  density  and  distribution  functions 
for  the  Rayleigh  distribution  are  defined  by: 


a.  Probability  density. 


^  lx  1  _  x  -X®  /2o® 

p  (  I  p  I  )  -  p  e  p 


(60) 


b.  Probability  distribution. 
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RATIO  = 


0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5 

INSTANTANEOUS  VALUE  PEAK  VALUE 

rms  VALUE  ‘  °  rms  VALUE 


(A)  Instantaneous  values  of  broad  band 
and  narrow  band  random  vibration  - 
Gaussian  distribution  (Equation  54); 

(B)  Peaks  of  broad  band  random  vibration 
(Equation  54); 

(C)  Peaks  of  narrow  band  random  vibration 
(Equation  61)  ; 

(D)  Instantaneous  values  of  a  sine  wave 
(Equation  55)  ; 

(E)  Peak  values  of  a  sine  wave. 


Figure  19.  Probability  Distribution  Functions. 
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(61) 


Since  the  Rayleigh  distribution  has  only  positive  values,  the  density 
function  can  also  be  written  as: 


P  cT 


(62) 


The  relations  given  by  equations  (60)  and  (61)  are  shown  graphically  by 
curve  A  of  Figure  20  and  curve  C  of  Figure  19  respectively.  Figures  18 
through  20  are  plotted  on  the  basis  that  the  mean  value  of  the  function  x(t) 
is  zero.  In  the  more  general  case  of  nonzero  mean  value,  the  absicssae  of 

these  figures  are  ^  ^  ^  ^  where  x  is  the  mean  value  of  the  function  x(t) . 

However,  in  vibration  data  analysis,  it  is  common  practice  to  constrain  the 
mean  value  to  zero  by  use  of  a  high-pass  filter  prior  to  statistical  analysis 
(elimination  of  d-c  mean).  Thus,  these  figures  show  the  form  usually 
encountered.  Figures  14  through  20  show  that  the  probability  density  and 
probability  distribution  curves  for  sinusoidal  and  random  functions  differ 
considerably,  thus  providing  identifying  characteristics  for-  each  type  of 
function.  When  the  time-history  x(t)  is  a  combination  of  sinusoidal  and 
random  functions,  the  shape  of  the  probability  density  and  distribution 
curves  depends  on  the  relative  magnitudes  of  each  type  of  function. 


Figure  20.  Normalized  Probability  Density  Functions. 
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2. 3. 3. 6  Gaussian  Distribution 


The  Gaussian  or  normal  distribution  mentioned  earlier  is  described 
by  the  density  function: 

-(x  -a)=/2o^ 
e 

=  on/^  (63: 

When  plotted,  it  has  the  characteristic  bell-shaped  curve  of  Figure  21. 

A  P(x) 


Figure  21.  Density  Function,  Gaussian  Distribution. 


The  curve  is  symmetrical  about  the  point  where  x  equals  a  and  has  a  width 
proportional  to  a.  It  can  be  shown  that  the  mean  value  of  the  distribution 
above  is  truly  a,  by  solving  for  the  first  moment.  Thus: 


=  J  X  p(x)  dx 
-00 

-  a)='/2o* 

=  - dx 


To  solve  this  integral,  a  change  of  variables  is  required  by  letting: 


y  =  ;  dx  =  V2o  dy 

oVT' 


then: 
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then : 


(V2„y  +  a)--y°,., 


y  e  dy  =  0 


J  e  ^dy  =  ''/tT 


MTP  5-1-025 
10  June  1968 


In  a  similar  manner,  it  can  be  shown  that  the  second  central  moment  or  variance 
ijs  is : 


Using  the  same  change  of  limits  as  described  in  Reference  23,  then; 


VT  a^y®  -y* 
- i—  e  dy 


4o*  r  2  -y®  , 


=  0^ 
Vtt  4 


Likewise,  it  can  be  shown  that  the  Gaussian  distribution  is  properly 
normalized  or  that; 
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J  p(x)  dx  =  1 

—00 

(74) 

oV2n  VtT 

-00 

(75) 

VTT 

O 

(76) 

This  Gaussian  curve  weighs  values  of  x  near  a  most  heavily.  The 

1 


value  of  p(x)  at  the  peak  is 


ov'2rr 


so  that  as  the  width,  cr,  decreases,  the 


height  of  the  curve  in  the  vicinity  x  =  a,  increases.  Ultimately,  as  o 
approaches  zero,  this  curve  approaches  the  delta  function  (x-a)  and  the 
variable  x  becomes  a  constant  a,  with  a  probability  of  1 .  It  is  often  de¬ 
sirable  to  know  the  relationship  between  the  true  rms  value  of  a  voltage 
which  is  normally  distributed  and  the  rectified  average  value  of  the  same 
voltage.  From  the  definition  of  the  second  moment,  it  has  been  shown  that 
the  true  rms  value  of  the  Gaussian  distribution  is  the  standard  deviation 
or  simply  a.  From  the  symmetry  of  the  distribution,  the  average  value  can 
be  found  by 


X 

av  ~ 


dx 


Substituting  y  =  — .  yields 

V2^ 


o 


=yi 


^  a  =  0.798  o  = 


(77) 


(78) 


(79) 


(80) 


Thus  the  ratio  of  rms  to  average  (form  factor)  of  a  normally  distributed 
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voltage  is  =  1 


=  1.25. 


It  has  already  been  shown  that  the  cumulative  probability  distri¬ 
bution  for  a  Gaussian  density  is  given  by; 


.(X)  -  J 


,-(x  -  a)®/2a® 


a/2r\ 


dx 


(81) 


and  can  be  expressed  in  terms  of  the  error  function  as: 


P(x)  -  h  y  erf  ( 


(82) 


The  average  value  of  a  Gaussian  function  cannot  be  defined  in  terms  of  a 
peak  value  since  the  peak  values  vary  from  zero  to  infinity.  However, 
physical  processes  are  usually  only  approximately  Gaussian  in  that  the  peak 
values  are  limited  to  some  value  between  three  and  ten  times  the  rms  value. 

2. 3. 3. 7  Rayleigh  Distribution 

It  has  been  stated  that  a  Rayleigh  distribution  arises  in 
vibration  work  when  the  peak  values  of  a  narrowly  filtered,  normally  dis¬ 
tributed  voltage  are  examined.  This  distribution  is  defined  by: 


P(x) 


-x=*/2cy^ 

xe _ 


(83) 


Here,  the  random  variable  x  refers  to  a  peak  value  of  the  filtered  noise 
signal  and  a  represents  the  rms  value  or  standard  deviation  of  the  original 
unfiltered  noise  signal.  Since  the  distribution  is  defined  for  positive  x 
only,  the  mean  value  or  first  moment  of  the  distribution  is  found  to  be  non¬ 
zero  and  is 

OO  CO 

m^  =  J  X  p(x)  dx  =  ~  J  ^  ^ 

o  o 


Then,  from  entry  509  of  Reference  23: 
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i.e.,  the  average  value  of  the  envelope  is  1-25  times  the  rms  of  the  original 
noise  wave.  Similarly,  the  second  moment  can  be  evaluated  as  follows: 

00 

ITlg  =  J  p(x)  dx 
o 

00 

_  1  r  3  -x®/2c^ 

-  J  X  e  dx 

o 


(86) 

(87) 


If  we  substitute  y  =  x^/2c^  then,  from  entry  508  of  Reference  23: 


00 

mg  =  2c^  J  y  e  ^  dy  =  2c^  (88) 

o 


That  is,  the  mean  square  of  the  envelope  is  twice  the  mean  square  of  the 
original  noise  wave.  Also  by  direct  substitution: 


^  (89) 

=  J  2^  -  =  V  0.429r*  (90) 

Thus  =  0.6557  or  the  rms  alternating  current  (a-c)  component  of  the 
envelope  is  0.655  of  the  rms  of  the  original  noise  wave. 

It  can  be  shown  that  the  Rayleigh  distribution  is  properly 
normalized  by  evaluation  of: 


X 


r 

J 

o 


p(x)  dx 


This  becomes: 


I 


xe 


-^I2d^ 


dx 


00 

p(x)  dx 
o 


which  reduces  to: 


00 


O 


then,  from  entry  429  of  Reference  23: 


an 


O 


The  probability  that  the  envelope  exceeds  some  specified  value 
by  integrating  the  density  function  from  x^  to  oo.  Thus: 


Prob. 


[envelope  >  “  J  P(x)  dx 

Xn 


This  is  plotted  in  Figure  22. 


Figure  22.  Probability  That  Envelope  Exceeds  x  . 

(Where  x  is  Multiple  of  a) 
o 
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is  obtained 
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It  is  interesting  to  note  that  the  Rayleigh  distribution  is 
independent  of  the  original  center  frequency  (about  which  narrow  band 
filtering  was  performed)  and  is  a  function  only  of  the  mean  square  value 
of  the  original  unfiltered  noise  wave.  Peak  factors  of  the  envelope  density 
function  have  been  partially  tabulated.  The  tabulation  is  listed  in  Table  1. 


Table  1.  Density  Function  Peak  Factors. 


%  Time  Exceeded 

Ratio 

Ratio  (db) 

10 

1.517 

3.62 

1 

2.146 

6.63 

0.1 

2.558 

8.39 

0.01 

3.034 

9.64 

0.001 

3.392 

10.61 

0.0001 

3.675 

11.31 

Here  the  peak  factor  ratio  is  given  as:  iL- 

This  table  is  for  the  peak  factors  of  the  Rayleigh  distribution  and  should 
not  be  confused  with  the  peak  factor  of  the  Gaussian  distribution  which  ap¬ 
plies  to  the  original  signal. 

2. 3. 3. 8  Determination  of  Distribution  Functions 

In  practice,  the  characteristics  of  a  signal  x(t),  at  a  specific 
time  and  at  a  specific  frequency  cannot  be  determined  since  only  the  in¬ 
stantaneous  value  of  x(t)  is  known  at  a  specific  time.  Instead,  the 
characteristics  of  the  signal  in  a  given  frequency  bandwidth  during  a  short 
period  of  time  are  determined.  Thus,  the  first  step  in  a  statistical  analysis 
is  to  select  a  time  interval  short  enough  so  that  the  nature  and  magnitude  of 
x(t)  may  be  assumed  constant,  but  long  enough  to  give  a  statistically  signifi¬ 
cant  result.  Often  it  is  convenient  to  store  the  data  sample  on  a  continuous 
loop  of  magnetic  tape,  so  that  it  may  be  played  back  repeatedly.  Then  a  time 
sample  of  the  signal  may  be  analyzed  to  determine  the  probability  distribution 
of  either  the  instantaneous  or  the  peak  values  of  the  signal.  The  concepts 
involved  in  this  determination  will  be  considered  separately.  The  instantaneous 
values  of  the  filtered  signal  are  indicated  by  x'(t)  while  the  peak  values  of 
the  filtered  signal  are  indicated  by  x^p(t).  The  instantaneous  values  and  the 
peak  values  are  described  as  follows: 

a.  Instantaneous  Values  -  The  probability  distribution  of 
instantaneous  values  is  obtained  for  a  set  of  discrete  values  of  x'(t) 
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rather  than  as  the  continuous  function  that  is  illustrated  in  Figure  19.  The 
probability  distribution  function  at  a  value  x is  the  probability  p(x">) 

that  X  (t)  exceeds  x  ,  .  (The  filtered  signal  x^(t)  may  be  passed  first  through 
a  mechanism  for  changing  the  polarity  of  the  signal  so  that  the  positive  and 
negative  values  of  x(t)  can  be  analyzed  separately  thus  halving  the  number 
of  discriminators  and  counters  required  for  a  given  resolution.  If  this  is 
done,  the  time -history  must  be  played  back  for  twice  as  long  but  this  usually 
causes  only  a  small  increase  in  the  to^tal  analysis  time.)  Refer  to  Figure  23 
and  note  that  this  is  equivalent  to  summing  all  the  At  ^  values  during  which 
X  (t)  exceeds  x  ,  or  level  L,  ,  and  dividing  by  the  total  elapsed  time.  If 

this  is  done  for  a  number  of  selected  levels,  Lj  ,  Ig  ,  ...  1^,,  the  probability 

distribution  function  at  these  levels  is  obtained.  Figure  24  illustrates  the 
block  diagram  of  a  statistical  analysis  measurements  system.  The  filtered 
signal  is  fed  to  an  array  of  discriminators  which  determine  whether  a  signal 
is  greater  or  less  than  a  preselected  value.  The  levels  of  the  discriminators 
are  set  at  ,  Ig  ,  L3 ,  etc.  usually  in  equal  increments  of  voltage.  The 

detail  or  resolution  obtained  is  governed  by  the  number  of  discriminator 
levels  selected  for  measurement.  The  detail  and  accuracy  yielded  by  ten 
levels  for  each  polarity  are  sufficient  for  many  engineering  purposes.  The 
resolution  obtained  improves  with  the  number  of  counters  (shown  in  Figure  24) 
registering.  Therefore,  a  variable  gain  control  is  needed  to  adjust  the 
signal  level,  by  a  known  amount,  until  an  adequate  number  of  counters  are 
registering.  When  the  signal  x^(t)  exceeds  the  level  of  a  particular  dis¬ 
criminator,  a  clock  is  started.  It  is  stopped  when  x^(t)  next  falls  below 
that  level,  and  the  counter  counts  the  elapsed  time  between  these  events. 

At  the  end  of  the  data  sample,  the  readings  of  the  counters  are  recorded 
by  the  read-out  mechanism  such  as  a  digital  counter.  If  each  reading  is 
divided  by  the  total  elapsed  time  of  the  data  sample,  the  values  obtained  are 
the  values  of  the  probability  distribution  function  at  the  values  of  the 
levels  Lq,  Lj ,  Lg ,  etc.  If  the  values  at  adjacent  levels  are  subtracted, 

the  average  probability  density  function  between  these  two  levels  is  obtained. 
The  polarity  of  the  signal  is  now  changed  and  the  process  is  repeated. 

However,  using  positive  polarity,  the  distribution  function,  p(x'>)  is 
obtained,  while  using  the  negative  polarity  yields  p(x^  <) .  Since 
P(x  <)  =  l-p(x  >) ,  the  two  sets  of  values  may  be  combined  to  give  the 
complete  distribution  function.  The  character  of  the  signal  can  then  be 
evaluated  by  the  use  of  probability  graph  paper. 

In  addition  to  determining  the  characteristics  of  the  signal,  this 
analysis  can  be  used  to  determine  the  mean  square  values  of  x^(t),  either 
by  use  of  the  probability  paper,  if  the  signal  is  Gaussian,  or  in  any  case, 
by  computing  the  moment  of  inertia  (second  moment)  of  the  probability  density 
function  about  the  mean  value  as  follows: 
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Figure  24.  Block  Diagrams  Showing  Steps  in  Statistical  Analysis. 
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*+•00 


0^  =  J  (x  -  x)^  p(x)  dx 
-00 


(96) 


=  ma  -  mi®  (97) 

b.  Peak  Values  -  The  probability  distribution  of  peak  values  or 
maxima  x"p(t)  of  the  filtered  signal  x"(t)  is  obtained  for  a  set  of  discrete 
values  of  x"p(t)  rather  than  as  the  continuous  function  illustrated  in 
Figure  19.  The  probability  distribution  function  at  a  value  x ^P| ,  is  the 

probability  p(x^Pi  s)  that  a  peak  or  maximum  of  x^(t)  exceeds  x^p^  .  Refer 

to  Figure  23  and  note  that  this  is  equivalent  to  counting  the  number  of  times 
that  x'’(t)  exceeds  x  ^p^  ,  (i.e.,  level  )  and  dividing  by  the  total  number 
of  times  that  x  (t)  exceeds  Lq .  If  this  count  is  made  for  a  number  of  levels 

Li  ,  Ig  ...  Lj^,  the  probability  distribution  function  at  these  values  is 

obtained . 


A  block  diagram  of  a  system  used  for  obtaining  signal  peak  measure¬ 
ments  is  shown  in  B  of  Figure  24.  The  filtered  signal  x^(t)  passes  through  a 
bank  of  level  discriminators  and  causes  the  associated  counters  to  register 
one  count  each  time  the  level  of  the  particular  discriminator  is  exceeded. 
Sometimes  an  array  of  discriminator-counter  combinations  is  called  a  "pulse- 
height  analyzer")  .  At  the  end  of  the  data  sample,  the  readings  of  the 
counters  are  recorded.  Division  of  each  count  by  the  zero  level  count 
yields  the  probability  distribution  function  p(x  p  S)  at  the  values 
equivalent  to  the  levels  ,  Lg  ,  ...  Ljj.  If  the  values  at  adjacent  levels 

are  subtracted,  the  average  probability  density  function  between  these  two 
levels  is  obtained.  The  character  of  the  signal  can  then  be  evaluated  by  the 
use  of  probability  graph  paper. 

For  some  end  uses  of  the  data,  such  as  fatigue  analysis,  the 
probability  distribution  of  the  total  or  unfiltered  signal  x(t)  rather  than 
the  probability  distribution  of  the  filtered  signal  x^(t)  is  desired.  Then, 
the  distribution  of  the  peaks  of  x(t)  is  obtained  by  shunting  out  the  filter 
as  indicated  by  the  dotted  lines  in  Figure  24. 

The  systems  shown  in  A  and  B  of  Figure  24  are  relatively  complex, 
and  unless  a  large  amount  of  analyses  of  this  type  are  to  be  made,  the  cost 
of  such  systems  may  not  be  justified.  A  less  complex  method  is  shown  in  C 
of  Figure  24.  In  this  case,  a  direct  writing  recorder  is  employed  and  the 
number  of  peaks  which  exceed  certain  levels  (such  as  shown  in  Figure  23)  are 
counted  manually.  In  principle,  the  distribution  of  instantaneous  values 
may  be  measured  in  this  way  but  practical  difficulties  make  this  method 
difficult  to  employ. 
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Recently,  techniques  have  been  developed  which  permit  estimation  of 
the  probability  density  or  distribution  function  of  broadband  random  signals 
with  accuracies  of  better  than  one  percent.  Signals  in  the  frequency  range 
of  d-c  to  10  kc  are  analyzed  by  means  of  a  slicer,  or  window,  circuit  which 
passes  pulses  of  a  radio  frequency  (RF)  oscillator  output,  normally  10  mega¬ 
cycles  (me)  every  time  the  signal  under  study  passes  through  the  specified 
window.  By  automatically  sweeping  the  window  levels  and  counting  the  cycles 
of  RF  passed  through  the  selecting  circuit  as  compared  to  the  total  number  of 
cycles  of  RF  passed  through  the  selecting  circuit  as  compared  to  the  total 
number  of  cycles  of  RF  generated,  accurate,  automatic  displays  of  probability 
density  within  the  range  of  ±  5  sigma  have  been  obtained. 

2.3.4  Frequency  Analysis 

The  objective  of  spectral  analysis  is  to  determine  the  variation 
of  vibration  magnitude  with  frequency.  The  narrowness  of  the  bandwidth  of 
the  filter  employed  in  the  analysis  determines  the  frequency  resolution  of  the 
analysis  and,_  therefore,  the  ability  to  detect  the  "fine-grain"  variation  of 
magnitude  with  frequency.  The  magnitude  obtained  at  a  particular  frequency 
will  be  the  average  magnitude  over  the  short  time  interval  of  data  analyzed 
during  which  the  nature  and  magnitude  of  the  vibration  may  be  considered 
constant.  (See  Figures  2  and  4.) 

In  preparing  for  spectral  analysis,  the  characteristics  of  the  time- 
history  x(t)  should  be  known  either  from  statistical  analysis  or  from  previous 
experience  so  that  the  most  appropriate  units  can  be  selected  for  describing 
the  vibration  magnitude.  In  practice,  the  spectral  analysis  often  is  con¬ 
ducted  first  and  the  necessity  for  statistical  analysis  is  determined  from  the 
resulting  spectrum. 

2. 3. 4.1  What  Value  to  Plot 

To  achieve  the  greatest  flexibility  from  vibration  data  and  make  it 
amenable  to  comparison  with  data  analyzed  on  other  systems,  it  should  be 
plotted  on  a  per-cycle  basis.  In  other  words,  the  quantities  used  to  express 
the  severity  of  vibration  in  different  parts  of  the  spectrum  should  be  divided 
by  the  effective  filter  bandwidth.  Further,  any  single  number  used  to 
describe  the  vibration  within  a  given  frequency  band,  or  in  the  entire  spectrum, 
must  be  obtainable  by  an  rms  operation,  from  a  plot  of  the  spectrum. 

In  the  determination,  of  which  single  number  quantity  should  be 
used  to  express  the  vibration  in  a  specified  band,  three  considerations  favor 
the  rms  of  the  instantaneous  accelerations. 

a.  The  rms  value  by  definition  is  the  standard  deviation  of  the 
instantaneous  accelerations  about  zero  as  a  mean  and  lends  itself  conveniently 
to  statistical  analysis. 

b.  The  rms  values  for  two  or  more  random  collections,  even  though 
they  are  in  different  frequency  bands,  generally  may  be  combined  as  the  square 
root  of  the  sum  of  the  squares . 
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c.  The  ms  value  is  the  widely  accepted  standard  measure  of  a-c 
voltage  in  electrical  circuitry. 

Plotting  the  function,  W(f)  against  frequency,  it  should  be  possible  to 
calculate  the  ms  acceleration,  o,  between  any  two  frequencies,  fj^  and  fg  , 

by  a  simple  operation,  based  on  the  requirement  that  the  ms  value  for  an 
entire  frequency  band  be  the  square  root  of  the  sum  of  the  mean  square  values 
for  all  bands  within.  This  is  satisfied  if  W(f)  is  defined  as: 


W(f)  = 


Lim  crAf® 
Af 


(98) 


where  oAf®  is  the  mean  square  acceleration  in  the  rectangular  bandwidths  Af> 

It  follows  then  that  the  function  W(f),  called  the  mean  square 
acceleration  density,  in  g^  per  cycles  per  second  (cps),  is  the  most  con¬ 
venient  means  of  describing  the  frequency  distribution  of  a  continuous 
spectrum. 

A  value  which  can  be  used  to  plot  frequency  distribution  of  power 
is  "power  spectral  density".  Power,  which  is  the  rate  of  doing  work,  is 
proportional  to  the  square  of  the  amplitude  of  a  hamonic  vibration.  If  two 
frequencies  are  present  in  a  vibration,  the  power  is  proportional  to  the  sum 
of  the  squares  of  the  individual  amplitudes  associated  with  the  two  fre¬ 
quencies. 


A  random  vibration  can  be  considered  to  be  a  sum  of  a  large  number 
(tending  to  infinity)  of  hamonic  vibrations  of  appropriate  amplitudes  and 
phase.  The  total  power  is  again  the  sum  of  the  component  hamonic  vibrations. 
The  way  in  which  this  power  is  distributed  as  a  function  of  frequency,  is 
known  as  "power  spectral  density"  and  it  represents  the  power  per  unit  fre¬ 
quency  interval.  A  plot  of  this  quantity  indicates  the  frequency  distribution 
of  power . 

2. 3. 4. 2  Definition  of  Power  Spectral  Density 

Power  spectral  density  is  defined  as  the  limiting  value  of  the  mean 
square  response  of  an  ideal  bandpass  filter  to  x(t),  divided  by  the 

bandwidth  B  of  the  filter  as  the  bandwidth  of  the  filter  approaches  zero. 

(An  ideal  bandpass  filter  has  a  transmission  characteristic  which  is 
rectangular  in  shape  so  that  all  frequency  components  within  the  filter 
bandwidth  are  passed  with  unity  gain  and  zero  phase  distortion,  while  fre¬ 
quency  components  outside  the  bandwidth  are  completely  removed.  If  the 
transmission  characteristic  is  H  instead  of  unity  in  the  bandwidth  B,  the 
spectral  density  is  obtained  by  dividing  the  mean  square  response  by  B  x 
instead  of  B.) 
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An  alternative  definition  is  as  follows:  If  the  function  x(t)  is 
passed  through  an  ideal  low  pass  filter  with  cutoff  frequency  f ,  the  mean 

square  response  of  the  filter  \  x'  will  increase  or  decrease  as  f^  is 

increased  or  decreased,  i.e.,  more  or  less  of  the  function  will  be  passed  by 
the  filter  (assuming  f^,  is  varied  in  a  frequency  range  where  the  power  spectral 

density  is  nonzero).  (An  ideal  low  pass  filter  is  an  ideal  bandpass  filter, 
having  a  lower  cutoff  frequency  of  zero.) 

The  power  spectral  density  W(f)  is  the  rate  of  change  of  [  \ 

with  respect  to  f^,,  i.e.. 


W(fc)  =  df7  t(x')=l 


(99) 


Figure  25  illustrates  a  plot  of  power  spectral  density  as  a 
function  of  frequency,  obtained,  for  example,  from  spectral  analysis  of  a 
random  function. 

# 


Figure  25.  Typical  Power  Spectral  Density  Plot  of 
Broadband  Random  Function  x(t). 


The  mean  square  value,  or  variance,  of  the  frequency  content  of  x(t)  between 
the  frequencies  f^  and  f^  is: 


0= 


(f. 


^  f  ^  fb)  = 


df 


(100) 
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_  The  ms  value  is  a-  The  mean  square  value  of  the  complete  function 

[  X  is  given  by  equation  (100),  when  f  and  f,  are  zero  and  infinity 

a  b 

respectively;  it  is  equal  to  the  area  under  the  entire  spectrum.  In  the  case 
of  white  noise,  for  which  the  spectral  density  is  independent  of  frequency, 
i.e.,  W(f)  =  W,  equation  (100)  simplifies  to: 

o*  =  W  (fg  -  fi)  (101) 


where  fg  and  fj^  are  the  limiting  frequencies  of  the  noise. 

To  understand  the  physical  meaning  of  the  power  density  spectrum, 
consider  one  simple  method  by  which  it  may  be  theoretically  measured. 

Suppose  that  the  power  density  spectrum  of  the  random  voltage  wave  fj^  (t) 

shown  in  Figure  1  is  required.  For  simplicity,  assume  that  the  random  wave 
contains  no  periodic  or  d-c  components.  To  detemine  the  spectrum,  the  random 
wave  is  applied  to  an  ideal,  conventional  type,  low  pass  filter  with  a  l-ohih 
resistor  output.  The  filter  cutoff  frequency  is  allowed  to  be  continuously  ad¬ 
justable  over  the  important  range  of  frequencies  of  the  random  wave.  Since 
the  filter  is  ideal,  all  components  of  the  random  wave  below  the  cutoff  fre¬ 
quency  pass  the  network  attenuated  and  all  others  are  completely  eliminated 
by  the  network.  By  means  of  an  appropriate  power  measuring  device,  the 
power  consumed  by  the  resistor  is  determined  at  various  cutoff  angular 
frequencies  ,  cife  ,  tea  and  plotted  as  shown  in  Figure  26.  Clearly,  a 

point  on  this  curve  at  a  given  angular  frequency  w  indicates  power  con¬ 
sumption  in  the  frequency  range  between  zero  and  the  given  frequency.  If 
the  slope  of  this  curve  is  determined  graphically  as  a  function  of  cu,  the 
result  would  appear  as  the  dotted  curve  shown  in  Figure  27.  This  curve  of 
slope  versus  angular  frequency  is  a  spectrum  of  power  density  because  its 
integral  from  zero  to  a  given  value  of  id  is  the  power  due  to  the  aggregate 
of  components  within  the  range  of  frequencies  under  consideration. 


CUTOFF  ANGULAR  FREQUENCY  OF  UOW-PASS  FILTER 

Figure  26.  Curve  of  Power  Consumed  by  Resistor  Versus  Cutoff 
Angular  Frequency  of  Ideal  Low-Pass  Filter. 
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Figure  27.  Slope  of  Curve  of  Figure  26  and  Power 

Density  Spectrum  of  Random  Wave  Figure  1. 


In  other  words,  the  measured  curve  of  Figure  26  is  an  integrated  power 
spectrum  so  that  its  derivative  as  given  in  Figure  27  is  a  power  density 
spectrum.  Finally,  in  order  that  the  power  density  spectrum  thus  obtained  be 
in  agreement  with  its  theoretical  definition,  the  amplitude  of  the  curve  for 
positive  frequencies  is  halved  and  the  upper  half  is  drawn  for  negative  fre¬ 
quencies  of  corresponding  values.  The  result  is  the  curve  6it(u;)  of  Figure  27 
which  is  an  even  function  as  it  should  be.  Negative  frequencies,  although 
physically  nonexistent,  have  been  introduced  in  harmonic  analysis  for  the 
reason  that  they  make  possible  the  expression  of  Fourier  series  and  integrals 
in  exponential  forms  which  are  simple  and  compact.  It  is  important  that 
expressions  involving  negative  frequencies  in  a  physical  problem  be  correctly 
interpreted.  Note  on  Figure  27,  that  where  '.u  equals  0,  the  curve  of 
has  a  finite  value;  however,  this  does  not  mean  that  the  random  voltage  has 
a  finite  d-c  component.  The  reason  is  that  inCiv)  equals  [W('x)J  is  a  spectrum 
of  power  density  so  that  the  d-c  power  according  to  equation  (116)  should  be 

W(o)  d,;' 

which  is  an  infinitesimal  quantity.  Actually,  this  amount  includes  the  con¬ 
tributions  of  all  components  in  the  band  d-.o  in  the  neighborhood  where  i 
equals  o.  Since  the  area  under  the  curve  of  ?,,(u;),  not  its  amplitude , 
represents  power,  the  amount  of  power  in  any  band  of  frequencies  finite  or 
infinite,  is  finite,  but  the  power  in  an  inf  initest  imal  band  is  in¬ 
finitesimal.  Furthermore,  the  power  at  a  particular  frequency  is  zero  when 
the  spectrum  is  a  continuous  density  spectrum. 

2. 3. 4. 3  Relationship  of  Power  Spectral  Density  to  Line  Spectrum 

The  mathematical  relationship  between  power  spectral  density  W(f) 

and  the  Fourier  coefficients  C  is; 

n 
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W(f)  =  2tt  Lim. 
B-»o 


jx'l^  =  2tt 
B 


L  CO'")’] 


=  %  ^  Cn®  6  (f  “  fn)  [b^  =  1] 
n=l 


(102) 


(103) 


where  [  x''  is  the  mean  square  response  of  the  ideal  filter  with  bandwidth 

B  and  center  frequency  f ,  where  f  equals  ^  ;  6  (f  -  fn)  is  the  Dirac  delta 
function. 

The  power  spectral  density  of  a  sinusoidal  function  is  a  line 

spectrum  since  such  a  function  has  effectively  zero  bandwidth.  [  x^  will 
change  instantaneously  as  the  cutoff  frequency  f^  increases  through  the  fre¬ 
quency  of  the  sine  wave  giving  theoretically  an  infinite  value  for  W(f). 

2. 3. 4. 4  Continuous  Spectrum  Versus  Line  Spectrum 

Full  information  concerning  harmonic  amplitudes  and  phases  is  con¬ 
tained  in  the  function  F(n),  which  in  general  is  complex  and  is  called  the 
complex  spectrum  of  f(t).  Therefore: 

+0O 

f(t)  =  ^  F(n)e^''“^  ^  (104) 

n=-oo 

F(n)  -i  .  dt  (105) 

2 


Because  the  harmonic  order  n,  assumes  only  discrete  values,  the 
spectrum  is  a  line  spectrum. 

Inasmuch  as  equation  (104)  is  a  summation  of  sinusoids  in  ac¬ 
cordance  with  the  amplitude  and  phase  information  contained  in  the  complex 
spectrum,  F(n),  for  regaining  the  original  function,  it  is  a  synthesis  in 
opposition  to  equation  (105)  which  is  an  analysis.  Again,  whereas  equation 
(105)  is  a  transformation  of  f(t)  into  its  frequency -domain  representation 
F(n),  equation  (104)  is  an  inverse  transformation  of  F(n)  into  its  time- 
domain  representation  f(t)  and  is  called  the  Fourier  transform  of  F(n). 
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Relationships  (104)  and  (105)  for  a  periodic  function  are,  therefore,  Fourier 
transforms  of  each  other  and  from  equation  (105)  we  get  a  complex  line 
spectrum.  If  equations  (104)  and  (105)  are  combined,  the  result  is: 


f(t)  Y.  do 

n=-03  ^ 

~2 

where  a  is  a  dummy  variable  introduced  for  the  sake  of  nonambiguity.  Now, 
if  the  period  grows  without  limit,  the  periodic  function  f(t)  tends  to  an 

aperiodic  one  and  in  so  doing,  equation  (106)  approaches  a  limiting  form.  As 
Tj  tends  to  infinity,  the  fundamental  angular  frequency  becomes  a 

differential  of  angular  frequency  dm  and  nOi  which  is  the  nth  harmonic  angular 

frequency,  becomes  the  continuous  angular  frequency  (.uj^  and  the  summation  over 

all  harmonics  becomes  an  integration  over  the  entire  continuous  frequency 
range  (-»,  <») .  Therefore,  the  limiting  form  for  equation  (106)  for  an 
aperiodic  function  f(t)  is: 


+03 


f(t)  = 


-oo 


-ho 


dm 


J 


da 


(107) 


This  is  the  Fourier  integral  for  an  aperiodic  function.  Writing 
two  reciprocals  from  equation  (107)  as  was  done  for  the  case  of  Fourier  series 
yields : 


f(t) 


+to 


r  F(m)e^^*’  dm 


F(m) 


-ho 


J  f(t)e'j''^ 


dt 


(108) 


(109) 


Here  F(m)  is  a  continuous  function  of  the  angular  frequency  m  and  is  in 
general  complex.  It  is  the  complex  continuous  spectrum  of  the  aperiodic 
function  and  the  relations  (108)  and  (109)  are  known  as  Fourier  transforms. 
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2.3.5  Determination  of  Power  Spectral  Density 

The  analysis  of  random  vibration  always  involves  a  record  of  finite 
length.  When  a  sample  of  time-history  of  finite  duration  is  employed  to  com¬ 
pute  the  power  spectral  density  of  a  random  function,  it  is  assumed  that: 

a.  The  function  is  ergodic,  i.e.,  that  averaging  one  time-history 
with  respect  to  time  yields  the  same  result  as  averaging  over  an  ensemble  of 
time-histories  at  a  given  instant  of  time. 

b.  That  the  function  is  stationary,  i.e.,  that  the  power  spectral 
density  is  independent  of  the  sample  of  time -history  chosen. 

Further,  the  averaging  time  or  sample  duration  must  be  long  enough 
to  yield  a  statistically  significant  value.  Thus,  the  mean  square  obtained 
should  not  vary  appreciably  with  a  change  in  averaging  time.  The  time  over 
which  a  vibration  record  may  be  considered  a  stationary  process  and  the  need 
for  a  sufficiently  long  averaging  time  often  are  conflicting  requirements. 

Assimie  that  a  magnetic  tape  record  of  finite  length  is  to  be 
analyzed  by  an  electronic  analyzer  as  a  stationary  process.  The  tape  can  then 
be  formed  into  a  loop  of  arbitrary  length  and  made  to  repeat  indefinitely  on 
the  analyzer.  By  such  a  procedure,  an  arbitrary  fundamental  period  correspond¬ 
ing  to  the  length  of  the  loop  is  established  and  the  contents  of  the  record  may 
be  defined  in  terms  of  the  Fourier  coefficients  of  integral  multiple  harmonics 
of  the  loop  frequency 

The  function  x(t)  then  can  be  represented  by  the  real  part  of  the 
Fourier  series  as: 


jnu/^t 


=  +  2  2^  jC^j  cos  (ncu^t  -  o/n)  (110) 

n=l 

where  the  complex  amplitude  is  defined  by  the  equation: 

T 

<^0  -  fi  ■  J  "y  (111) 

-T 

and  2T  is  the  loop  period. 

The  mean  square  value  x®  is  a  stationary  property  of  x(t)  that 
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is  determined  by  equation  (6)  as: 


x= 


T 

[x(t)?  dt 

t 

-T 


n=l 

This  equation  indicates  that  the  contribution  to  the  mean  square 
value  in  any  frequency  interval  is  the  summation  of  all  the  components  within 
that  frequency  band . 

If  now  the  length  of  the  loop  and  the  loop  period  are  doubled,  the 
number  of  spectral  lines  will  be  doubled.  Since  the  mean  square  value  of  the 
entire  spectrum  must  remain  essentially  constant,  the  values  of  the  new 

coefficients  must  decrease  to  approximately  one  half  their  former  values 

as  illustrated  in  Figure  28. 


w  /2  (B)  FREQUENCY 

o  _ 


Figure  28.  Spectrum  Analysis  of  Random  Vibration, 

As  Determined  from  a  Magnetic  Tape  Loop. 
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The  sum  of  the  lengths  of  the  lines  up  to  any  frequency  is  the  same  in 

each  diagram  of  Figure  28,  and  the  increment  in  x®  divided  by  the  increment 
in  frequency  is  essentially  the  same  in  each  case.  Thus: 


(113) 


The  quantity  is  the  power  spectral  density  at  u)  =  UJfj'  "^he 

dimensions  are  the  square  of  the  parameter  represented  by  in  equation  (113) 
per  unit  of  frequency. 

As  the  length  of  the  loop  is  increased,  more  and  more  spectral 
lines  are  introduced;  the  magnitude  of  jCul*  decreases  correspondingly  but 

W(G:f,)  remains  finite.  In  the  limiting  case,  as  T  approaches  infinity  or  Act 

approaches  zero,  the  spectrum  becomes  continuous  and  the  discrete  values  of 
W(ct^)  approach  a  smooth  function  W(ai): 


Ax^  dx* 
u  =  lim  — —  =  — — 

Att)  dcu 


(114) 


The  mean  square  value  can  be  expressed  by  the  integral: 


=  J  W(lo) 


(115) 


where  only  positive  frequencies  are  considered.  This  is  equivalent  to  the 
definition  over  the  entire  range  of  frequencies  given  by 


lim  1 
T-teo  2 


“  J  [fi(t)]^  "  J 


(116) 


The  power  spectral  density  of  a  random  time  function  may  be  deter¬ 
mined  by  feeding  the  function  x(t)  into  a  spectrum  analyzer  that  transmits 
only  those  frequency  components  within  the  passband  of  the  analyzer. 
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u;  ±  as  shown  in  Figure  29.  The  output  in  the  passband,  namely,  is 

2 

then  indicated  by  a  "mean  square"  meter.  The  mean  power  spectral  density 
W(u))  over  the  passband  at  tB  is  determined  by  dividing  the  meter  reading  by 
Au).  This  division  may  be  incorporated  directly  in  the  calibration  of  the 
meter.  The  complete  distribution  of  W((ju)  then  is  determined  by  changing 
the  frequency  u)  in  increments  of  The  mean  square  and  averaging  methods 

often  used  to  determine  power  spectral  density  are  shown  in  Figures  32  and  33. 

2. 3. 5.1  Measurement  of  Power  Spectral  Density  by  the  Use  of  Filters 

When  a  random  function  is  applied  to  the  input  of  a  filter,  the  out¬ 
put  is  not,  in  general,  represented  by  the  normally  defined  3-decibel  (db) 
bandwidth.  Rather,  the  "noise  bandwidth"  of  the  filter  is  involved  and  the 
output  is  obtained  by  multiplying  the  input  spectral  density  by  the  square  of 
the  transmission  characteristic  of  the  filter.  An  ideal  bandpass  filter  has 
a  bandwidth  B  and  an  output  spectral  density  of  W  times  B  if  W  is  the  constant 
spectral  density  input.  In  a  typical  case  however,  the  effective  or  noise 
bandwidth  of  a  filter  should  be  determined  by  the  mean  square  response  of  the 
filter  to  known  white  noise.  If  a  random  signal  is  applied  to  an  ideal  filter, 
the  power  spectral  density  output  is  the  mean  square  output  divided  by  B .  In 
a  practical  case,  the  mean  square  output  divided  by  the  effective  filter 
bandwidth  is  a  good  approximation  of  the  average  spectral  density  within  the 
filter  bandwidth. 


Figure  29.  The  Experimental  Determination  of  The  Power 

Spectral  Density  W(u;)  of  a  Random  Function  x(t). 
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2. 3. 5. 2  Effects  of  Filter  Characteristics 

The  properties  of  filters  used  for  the  measurement  of  power 
spectral  density,  such  as  transmission  pattern  characteristics,  time  constant, 
and  the  rate  of  variation  of  spectral  density  with  frequency  of  the  signal 
under  study  will  have  a  most  significant  effect  on  the  resulting  spectrum 
plot.  If  the  filter  bandwidth  is  wider  than  the  most  narrow  resonance  under 
study,  then  although  analysis  time  might  be  reduced,  there  is  insufficient 
resolution  for  an  accurate  analysis.  In  general,  the  effect  of  too  wide  a 
filter  bandwith  is  to  smooth  the  resulting  plot.  The  output  may  in  some 
cases,  be  time  averaged  by  the  effect  of  the  filter  time  constant.  The 
effect  of  too  wide  a  filter  (dotted  lines)  and  how  it  can  mksk  peaks  and 
notches  is  shown  in  Figure  30. 


Figure  30.  Effects  of  Filter  Bandwidth  Characteristics  on 
Measured  Spectra, 

2. 3. 5. 3  Measurement  of  Line  Spectra  by  the  Use  of  Filters 

Sinusoidal  functions  are  analyzed  with  filters  by  plotting  filter 
output  amplitudes  instead  of  spectral  density.  The  accuracy  of  the  measure¬ 
ment  depends  critically  on  the  filter  bandwidth,  and  on  the  readout  device. 
Two  Sinusoids  narrowly  separated  in  the  spectrum  are  detectable  only  with 
sufficiently  narrow  filtering  devices  (high  selectivity).  For  a  random  wave 
applied  to  a  narrow  band  filter,  the  output  is  often  considered  to  be  the 
mean  square  value  of  a  single  sinusoid  rather  than  the  mean  square  value  of 
the  sum  of  many  sinusoids. 

2. 3. 5. 4  Spectrum  of  Sinusoid  Plus  Noise 

’  ■  if  the  random  wave  under  observation  has  a  periodic  component,  then 
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either  the  power  spectrum  or  power  density  spectrum  could  prove  inadequate. 
Strictly  speaking,  the  power  density  spectrum,  although  well  suited  for 
characterizing  a  random  wave,  does  not  give  a  convenient  representation 
of  a  periodic  component  because  of  the  infinite  discontinuities  involved 
at  the  "line  spectrum"  points  of  the  periodicity.  Similarly,  the  random 
component  is  not  easily  represented  by  the  power  spectrum.  For  example,  in 
Figure  31,  the  peak  could  be  caused  either  by  a  periodicity  or  a  truly  random 
function  of  increased  intensity  over  that  frequency  band.  In  such  cases,  a 
periodicity,  if  present,  can  be  easily  detected  in  the  probability  density 
function  of  the  signal  under  study.  Presence  of  a  periodicity  can  also  be 
detected  in  the  integrated  power  spectrum,  defined  as: 

w 

J  W(co)  diu  (117) 

•<o 


This  function  increases  monotonically  for  a  random  wave  but  has  finite  dis¬ 
continuities  at  the  harmonic  frequencies  of  a  periodic  wave. 

2. 3. 5. 4  Measurement  of  Filter  Output 

The  magnitude  of  the  filter  output  can  be  obtained  in  terms  of  the 
mean  square  or  average  value  of  the  filtered  signal.  If  statistical  analysis 
has  shown  x'’(t),  the  filtered  signal,  to  be  a  random  function  then  the  mean 
square  value  is  determined  by  squaring  and  integrating  the  function  and  then 
dividing  by  integration  time  T  as  shovm  in  Figure  32.  Division  of  the  output 
of  the  integrator  or  alternate  averaging  device  (shown  dotted)  by  the  effective 
filter  bandwidth  yields  the  average  power  spectral  density  within  the  filter 
bandwidth.  When  the  filter  output  is  simply  rectified  and  smoothed,  the 
average  value  of  x'(t)  is  obtained.  If  this  is  known  to  be  a  random  function, 
power  spectral  density  can  be  obtained  by  the  relationships  described  in 
paragraph  2.3.3  or  by  using  the  empirical  measurement  of  response  to  a  known 
spectral  density.  Spectral  analysis  is  accomplished  by  the  averaging  method 
shown  in  Figure  33. 

2.3.6  Magnitude  Analysis 

2. 3. 6.1  Magnitude  Time-History  Analysis 

In  this  type  of  analysis,  the  variation  of  the  vibration  severity 
as  a  function  of  time  is  examined.  The  result  may  appear  as  in  Figure  34. 

Often  the  magnitude  is  computed  as  the  average  over  a  few  seconds  while  the 
data  record  itself  may  last  for  several  minutes  .  Either  the  magnitude  of  the 
complete,  unfiltered  signal  may  be  examined,  such  as  the  overall  mean  square 

value  X* ,  or  the  magnitude  of  the  vibration  in  a  restricted  band  of  frequencies 
can  be  determined.  If  a  spectrum  analysis  shows  that  the  important  components 
of  a  given  vibration  signal  exist  only  over  selective  frequency  intervals, 
then  the  magnitude  time-history  analysis  would  be  carried  out  only  for  those 
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FREQUENCY  f 


Figure  31.  Example  of  Line  Spectrum  Superimposed 
on  Random  Vibration  Spectrum. 


TYPICAL  WAVEFORMS 


WAVEFORM 


Figure  32.  Block  Diagram  Showing  Spectral 

Analysis  by  the  Mean-Square  Method. 
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frequencies.  For  example,  the  amount  of  analysis  necessary  to  define  the 
surface  of  Figure  2  depends  on  the  degree  of  the  completeness  with  which  the 
surface  is  to  be  defined. 


TYPICAL  WAVEFORMS 


x(t)  x'{t)  [*'] 


Figure  33.  Block  Diagram  Showing  Spectral 
Analysis  by  Averaging  Method. 


Figure  34.  Typical  Plot  of  Magnitude  Time- 
History  Analysis. 

2. 3. 6. 2  Magnitude  Distribution  Function 

The  magnitude  distribution  function  of  a  signal  is  the  total  time  that 
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the  magnitude  is  equal  to  or  greater  than  a  particular  magnitude.  The  mean 
square  or  mean  value  of  the  signal  is  measured  on  a  bank  of  discriminator- 
counter  or  discriminator -clock-counter  combinations,  dependent  on  the 
averaging  used.  (See  Figure  35.)  Figure  36  illustrates  one  method  of 
plotting  the  readings  from  the  counters.  The  ordinate  corresponds  to  the 
magnitude  of  the  level  which  is  exceeded  whereas  the  abscissa  is  normalized 
to  indicate  the  percentage  of  total  time  that  the  level  is  exceeded.  The 
function  illustrated  in  Figure  36  is  analogous  to  probability  density  as 
shown  in  Figure  19.  Unless  the  analysis  represents  the  total,  unfiltered 
signal,  the  filter  center  frequency  and  bandwidth  must  be  specified.  The 
relations  between  times  of  occurrences  of  different  vibration  magnitudes  is 
lost  with  this  type  of  presentation,  similar  to  the  way  time  is  averaged 
out  of  spectrum  presentation  on  all  but  real  time  analyzers  as  described  in 
Reference  30. 

2.3.7  Correlation  Analysis 

A  correlation  function  describes  the  time  relation  between  two 
parameters  as  a  function  of  the  time  delay  between  the  parameters  at  the 
times  in  which  they  are  observed.  Physically,  correlation  involves  the  delay 
of  one  signal  with  respect  to  another,  multiplication  of  the  two  signals  and 
integration  of  the  product  over  suitable  limits.  If  the  two  signals  represent 
the  same  phenomena,  but  delayed  in  one  case  with  respect  to  itself,  then  the 
process  is  called  autocorrelation.  For  example,  x^  (t)  and  x^  (t  +  t) ,  where 

XjL  (t)  is  a  physical  event  such  as  missile  vibration  at  a  specific  point.  If 

the  two  functions  are  physically  distinct  from  one  another  (for  example, 

Xj^  (t)  and  Xg  (t  +  t)  then  the  process  of  delay,  multiplication,  and  summation 

is  known  as  crosscorrelation. 

Correlation  analysis  can  be  applied  to  missile  vibrations  in  several 
ways.  For  example,  the  propagation  of  vibration  through  a  structure  can  be 
studied  by  examining  the  crosscorrelation  function  between  acceleration 
measurements  at  two  different  points  of  the  structure.  Crosscorrelation 
between  parameters  with  diffei'ent  physical  units  is  also  possible;  as  for 
example  the  crosscorrelation  between  applied  force  and  resulting  acceleration 
at  points  of  a  missile. 

2. 3. 7.1  Autocorrelation  Function 

The  autocorrelation  function  of  a  signal  x(t)  also  known  as  the 
average  lag  product  R(t)  is  defined  as: 

T 

lim  1  p 

— J  x(t)  x(t  +  t)  dt  (118) 

o 
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Figure  36.  Distribution  Function  of  Vibration  Magnitude 
for  a  Filter  Having  a  Center  Frequency  f 


and  is  seen  to  be  a  function  of  the  time  delay  t-  For  strictly  statistical 
usage,  division  by  the  normalizing  factor  R(o)  is  necessary,  so  that  the 
value  of  the  function  is  unity  at  zero  lag.  Under  this  condition,  it  can 
be  called  the  autocorrelation  function.  However,  widespread  use  has  brought 
about  the  name  autocorrelation  function  for  R(t)  and  not  the  ratio  R(t)/R(o), 


function : 


Several  useful  properties  ar4  associated  with  the  autocorrelation 


a.  R(t)  =  R(“t)  i.e.,  it  is  always  an  even  function. 


b.  R(o)  =  0*  +  i.e.,  it  equals  the  average  power  in  a  sample 

function  or  the  sum  of  a-c  and  d-c  powers. 

c.  R(o)  s  R(t)  i.e.,  the  largest  value  of  the  autocorrelation 
function  occurs  at  zero  lag  since  a  function  is  never  more  closely  related 
to  another  function  than  it  is  to  itself. 

d.  |R(t) I  S  1  i.e.,  the  autocorrelation  function  can  never  exceed 
plus  or  minus  one  in  magnitude. 

e.  The  Fourier  transform  of  the  autocorrelation  function  is  the 
power  density  spectrum.  In  other  words,  if  W(f)  is  the  one-sided  power 
spectral  density  then: 


R(t)  -  J  W(f)  cos  (JOT 


(119) 


R(t)  cos  u)t  dT 


(120) 
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The  last  two  equations  are  known  as  the  Wiener-Khint chine  relations 
and  form  a  Fourier  transform  pair. 

Equation  (120)  can  be  used  to  determine  the  power  spectral  density 
of  a  random  function  from  the  autocorrelation  function. 

f.  R(  oo)  =  0,  if  x(t)  contains  no  periodic  or  d-c  components. 

g.  If  the  autocorrelation  function  is  continuous  at  the  origin, 
it  is  continuous  everjrohere . 

One  of  the  most  important  characteristics  of  the  autocorrelation 
function  is  that  even  though  retaining  all  the  harmonics  of  the  given  function 
and  containing  no  new  ones,  it  discards  all  their  phase  angles.  In  other 
words,  all  periodic  functions  having  the  same  harmonic  amplitudes  but 
differing  in  their  initial  phase  angles,  have  the  same  autocorrelation  function. 
Thus,  the  power  spectrum  of  a  periodic  function  is  independent  of  the  phase 
angles  of  its  harmonics. 

2. 3. 7. 2  Determination  of  Autocorrelation  Function 

One  method  of  obtaining  a  correlation  analysis  by  analog  means  is 
illustrated  in  Figure  37.  An  almost  identical  arrangement  is  used  for  both 
autocorrelation  and  crosscorrelation  analysis.  For  autocorrelation,  a  stored 
signal  is  applied  to  a  multiplier  both  directly  and  through  a  variable  time 
delay  (t).  For  crosscorrelation,  two  separate  signals  are  applied  to  the 
multiplier,  one  delayed  t  seconds  with  respect  to  the  other.  A  plot  of 
normalized  integrator  output  versus  time  delay  is  the  autocorrelation  function. 

If  the  signal  is  available  in  digital  form,  a  standard  digital 
computer  can  be  programmed  to  rapidly  perform  correlation  analysis. 

2. 3. 7. 3  Crosscorrelation  Function 

The  crosscorelation  function  for  two  distinct  time  varying  parameters 
is  defined  as: 


T 

Ria  (t)  =  It  J  Xs  (t  +  t)  dt 

-T 


(121) 


and  is  again  a  function  of  t-  |Ria  (t) |  is  then  a  measure  of  the  correlation 

between  (t)  and  Xg  (t)  vdien  observed  at  any  two  times  separated  by  the  time 

difference  t-  Note  that  (t)  =  Rsi(-t)  and  that  R^g (l®)  =  0.  The  Fourier 

transforms  relating  the  crosscorrelation  function  to  the  cross -power  density 
spectrum  are  defined  as: 
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Figure  37.  Correlation  Analysis. 
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where  Wjs  (u))  is  the  cross -power  density  spectrum  of  the  random  functions 
Xj (t)  and  Xg  (t)  . 

2. 3. 7. 5  Determination  of  Crosscorrelation  Functions 

A  method  of  computing  crosscorrelation  functions  by  analog  means  is 
shown  in  B  of  Figure  37.  Digital  computation  techniques  can  be  used  to  good 
advantage  in  the  calculation  of  crosscorrelation.  In  either  method  the 
signal  Xj^  (t  +  t)  used  previously  for  autocorrelation  analysis,  is  replaced 

Xg  (t  +  t)  • 


and 

Ris(t)  = 
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3.0  CONCEPTS  IN  SHOCK  DATA  ANALYSIS 


In  connection  with  environmental  testing  or  structural  evaluations, 
shock  may  be  defined  as  a  pulse,  (transient  or  change  of  acceleration)  which 
stands  out  distinctly  above  the  general  vibration  level  and  terminates  in  a 
time  that  is  short  by  comparison  with  the  shortest  decay  time  associated  with 
the  structure  which  produced  the  data. 

3.1  Engineering  Uses  of  Shock  Measurements 

In  general,  a  shock  measurement  in  the  form  of  a  time-history  of 
the  acceleration  or  some  other  parameter  is  not  the  most  useful  type  of 
presentation  upon  which  to  base  an  engineering  evaluation.  The  measured 
shock  often  includes  the  vibrational  response  of  the  structural  element  to 
which  the  shock  measuring  device  is  attached.  This  vibration  could  obscure 
the  determination  of  the  actual  shock  phenomena  whose  presentation  is  desired. 
Data  reduction,  in  this  case,  could  eliminate  many  irrelevancies  of  the 
measured  data.  What  is  most  often  desired  is  a  presentation  from  which 
theoretically  predicted  values  can  be  compared  with  the  measured  results  or 
else  a  presentation  which  can  be  used  for  the  calculation  of  structural  re¬ 
sponse  . 

3.1.1  Concepts  of  Shock  Data  Reduction 

One  of  two  basic  methods  normally  is  employed  to  describe  shock. 

It  can  be  described  in  terms  of  the  inherent  properties  of  the  shock,  either 
in  the  time  domain  or  the  frequency  domain,  or  a  description  can  be  made  of 
the  shock  in  terms  of  its  effect  on  structures  when  the  shock  acts  as  the 
excitation.  The  applicable  method  will  be  determined  by  the  intended  usage 
of  the  reduced  data.  Regardless  of  the  type  of  reduction  employed,  a  record 
of  the  original  time-history  should  be  included  with  the  data  whenever 
possible . 


3.1.2  Reduction  to  the  Frequency  Domain 

Any  transient  function,  such  as  a  pulse,  can  be  broken  down  into  a 
continuous  Fourier  spectrum.  If  the  transient  is  recorded  on  a  tape  loop  and 
played  back  continuously,  a  good  approximation  of  the  Fourier  spectrum  can  be 
obtained  with  almost  any  wave  form  analyzer.  The  time  functions  and  resulting 
spectrum  functions  for  a  dozen  representative  transients  are  given  in  Figures 
38  through  40.  Graphical  analysis  for  the  determination  of  Fourier  coeffi¬ 
cients  can  often  be  employed  successfully,  however,  the  time  involved  in 
graphical  determination  as  well  as  the  availability  of  short  transient  auto¬ 
matic  recording  and  analysis  devices  and  digital  computers,  usually  results 
in  a  preference  for  electronic  analysis. 

The  description  of  a  signal  by  its  Fourier  spectrum  is  most  often 
associated  with  linear  system  analysis.  In  the  reduction  and  usage  of  shock 
data,  if  the  properties  of  the  structure  upon  which  the  shock  acts  are  defined 
as  a  function  of  frequency  (the  structure  is  often  assumed  linear  over  some 
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range  of  parameters)  they  are  generally  called  transfer  functions.  Commonly 
used  transfer  functions  in  shock  and  vibration  work  include  mechanical  im¬ 
pedance,  mobility,  and  transmissibility • 

The  transfer  function  of  a  structure  can  sometimes  be  determined  by 
applying  a  pulse  of  force  to  the  structure  and  observing  its  response.  If 
such  a  transfer  function  is  determined  to  be  T(u))  then  the  Fourier  spectrum 
of  the  response,  R(u)),  is  the  product  of  the  Fourier  spectrum  of  the  shock 
F((u)  and  the  transfer  function  T(u)). 

3.1.3  Reduction  to  the  Response  Domain 


Vfhen  a  shock  is  applied  to  a  structure,  a  response  characteristic 
of  the  structure  ordinarily  is  produced-  If  certain  assumptions  concerning 
the  properties  of  the  structure  are  made,  the  response  can  be  assumed  to 
depend  only  on  the  damping  of  the  structure.  The  magnitudes  of  the  peaks  of 
the  response  can  then  be  systematically  related  to  the  properties  of  the 
system.  The  description  of  a  shock  in  terms  of  its  Fourier  spectrum  differs 
from  this  in  two  ways : 

a.  If  a  knowledge  of  the  peak  responses  of  a  structure  to  a  shock 
excitation  exists,  the  time-history  of  the  shock  still  cannot  be  determined 
on  this  basis  alone;  i.e.,  the  calculation  of  the  peak  responses  is  an  ir¬ 
reversible  process. 

b.  The  response  spectrum  describes  the  effect  of  a  shock  upon  a 
structure  in  terms  of  its  peak  responses  where,  on  the  other  hand,  the  Fourier 
spectrum  defines  a  shock  in  terms  of  the  amplitudes  and  phases  of  its  fre¬ 
quency  components. 

If  a  representation  of  the  peak  responses  of  many  simple  structures 
is  obtained  from  the  shock  measurement,  the  process  is  called  reduction  to 
the  response  domain.  For  a  system  that  responds  in  a  single  degree  of  freedom 
this  type  of  reduction  applies  directly.  However,  for  a  system  that  responds 
in  more  than  one  degree  of  freedom,  this  method  of  reduction  is  applicable 
only  if  the  principle  of  superposition  can  be  used  with  validity. 

3.1.4  The  Shock  Spectrum 

The  shock  spectrum  for  a  given  shock  is  defined  as  the  peak  accelera 
tion  response  of  a  simple  undamped  resonant  member  to  this  shock,  plotted  as  a 
function  of  the  frequency  to  which  the  member  might  be  tuned.  It  thus  serves 
as  an  indication  of  possible  stress  within  the  component,  whereas  the  peak 
acceleration  of  the  original  complex  serves  merely  as  an  indication  of  stress 
at  the  mounting  points  of  the  component. 

In  common  terminology,  the  term  shock  spectrum  is  used  synonymously 
with  two-dimensional  shock  spectrum  and  represents  the  maximum  value  of  the 
response  found  in  a  single  time-history.  The  three-dimensional  shock  spectrum 
is  also  encountered  and  takes  the  form  of  a  surface  which  shows  the  distri¬ 
bution  of  response  peaks  throughout  the  time-history. 
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3. 1.4.1  Relation  Between  Shock  Spectrum  and  Fourier  Spectrum 
If  F(ya)  is  the  Fourier  spectrum  of  a  pulse  where: 


+00 


F((ju)  =  J  F(t)e"^“^^  dt 

mCD 


(124) 


then  for  a  pulse  starting  at  t  equals  0,  it  is: 


F(u))  =  J  f(t)e‘^‘“*'  dt  (125) 

o 

It  can  be  shown  that  a  quantitative  relationship  exists  between  the 
Fourier  spectrum  of  absolute  values  and  the  peak  residual  response  of  an  un¬ 
damped  simple  structure. 

If  F(oJj^)  is  the  Fourier  amplitude  spectrum  or  the  absolute  value  of 

F(m)  at  the  frequency  then  an  equivalent  static  acceleration  (A  eq.)  can 

be  assigned  to  the  maximum  residual  displacement  response.  If  the  fraction 
of  critical  damping  of  the  simple  structure  used  for  response  measurements  is 
zero  then : 


A  eq. 


■— n-i.) 


(126) 


where  g  is  the  acceleration  due  to  gravity. 

3.1.5  Methods  of  Data  Reduction 

Generally,  in  shock  data  reduction,  the  measured  shock  is  not  de¬ 
finable  by  an  analytic  function  but  is  available  in  graphic  or  electrical 
form.  Several  reduction  methods  have  been  described  in  preceding  paragraphs. 
However,  a  general  summation  of  available  methods  of  shock  data  reduction 
would  have  to  include  reduction  to  the  Fourier  spectrum  (by  wave  analysis, 
digital  computation,  use  of  trigonometric  forms,  and  tables)  and  reduction 
to  the  shock  spectrum  (by  graphic  methods,  numerical  methods,  mechanical, 
electrical  analog  methods,  and  use  of  the  reed  gauge). 

4.0  CONCEPTS  OF  ACOUSTICAL  DATA  ANALYSIS 

4 . 1  General 
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A  qualitative  measure  of  the  relative  "loudness"  of  noise,  or  some 
other  sound,  is  often  made  in  connection  with  the  testing  and  evaluation  of 
jet  aircraft  and  missile  systems.  This  measurement  can  then  be  used  for  the 
evaluation  of  personnel  safety,  system  performance,  or  specification  design. 
With  the  recent  introduction  of  large,  reverberant  test  chambers  and 
sophisticated  equalization  equipment,  acoustical  testing  can  conceivably  be 
performed  directly  from  the  raw  data,  much  in  the  same  way  that  inflight 
vibration  data  can  conceivably  be  faithfully  reproduced  on  a  shaker  table. 
Thus,  in  vibration  and  acoustics,  a  measure  of  the  amplitude  of  the  environ¬ 
ment  as  well  as  its  frequency  spectrum  will  serve  many  useful  purposes . 

4.1.1  Units  of  Measure 


In  the  description  of  the  loudness  of  sound  many,  sometimes  con¬ 
fusing,  terms  are  used  and  quite  often  coincidentally.  These  terms  include 
sones,  phons,  decibels,  spectrum  level,  and  the  relative  term  -  loudness. 

So  that  a  scale  of  loudness  as  a  function  of  the  sound  pressure 
level  in  db  can  be  obtained,  subjective  measurements  are  made  with  a  pure  tone. 
The  frequency  of  this  tone  is  normally  1000  cps  and  a  reference  level  is 

chosen  at  40  db  referenced  to  (re)  2  x  10  ^  U  bar  and  is  equal  to  one  sone  or 

1000  millisones.  The  intensity  of  sound  normally  is  expressed  in  terms  of 
sound  pressure  in  p,  bars  and  is  equal  to  one  dyne  per  centimeter  squared 

(dyne/cn^).  A  loudness  of  10  sones  is  then  ten  times  as  loud  as  1  sone.  The 

relationship  between  loudness  in  millisones  and  loudness  level  in  phons  has 

been  determined  empirically  and  produces  a  curve  which  shows  that  at  1000  cps, 
the  loudness  level  in  phone  corresponds  to  the  sound  pressure  level  in  db  re 

2  X  10  ^  IX  bar.  The  dependence  of  loudness  level  (phons)  on  frequency  has 
been  established.  Thus,  loudness  in  sones  can  be  applied  to  sound  of  any  fre¬ 
quency  or  combination  of  frequencies.  As  an  example  of  the  non linearities 
involved,  a  reduction  in  loudness  level  from  70  phons  (10,000  millisones)  to 
40  phons  (1000  millisones)  corresponds  to  a  tenfold  change  in  loudness.  But 
a  change  from  40  phons  to  20  phons  is  again  a  tenfold  change  in  loudness  to 
100  millisones . 

4.1.2  Addition  of  Loudness  Levels 


Because  of  phasing,  masking,  and  other  considerations,  loudness  in 
general  is  not  directly  additive.  Since,  however,  the  frequency  components 
of  a  given  noise  are  often  determined  by  filtering  the  overall  noise,  it  has 
been  found  convenient  to  express  total  loudness  of  data  passed  through  nth 
octave  bandpass  filters.  This  results  in  the  equation: 


=  S  +  F(i:s  -  S  ) 
t  m  m 


(127) 
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Where  S  is  the  total  loudness  in  sones,  S  is  the  loudness  in  sones  of  the 
t  m 

loudest  nth  octave  band,  and  2s  is  the  sum  of  the  loudness  in  sones  of  all  the 
nth  octave  bands .  By  taking  the  bandwidth  and  masking  effect  into  account, 
the  following  F  factors  were  obtained: 

a.  1/3  octave  bandwidth  where  F  =  0.15 

b.  1/2  octave  bandwidth  where  F  =  0.2 

c.  1  octave  bandwidth  where  F  =  0.3 

This  is  also  shown  graphically  in  Figure  41. 

In  addition  to  equal  loudness  contours,  the  concept  of  equal  noisi¬ 
ness  contours  has  also  received  much  attention  and  gives  a  comparable  answer 
in  the  calculation  of  loudness  level  of  a  specific  sound. 


1/3  1/2  I 

WIDTH  OF  BAND  IN  OCTAVES 


Figure  41.  Value  of  F  for  Bands  of  Various  Widths. 

4.1.3  Pure  Tones  -  Discrete  Frequencies 

A  pure  tone  in  acoustical  data  analysis  is  analogous  to  a  vibrational 
acceleration  at  a  discrete  frequency  and  both  result  in  a  line  spectrum. 

Sounds  with  line  spectra  can  be  analyzed  with  any  type  of  analyzer  which  has 
sufficient  resolution.  The  presentation  of  the  data  may  be  either  in  tabular 
form  (spectral  linesversus  frequency)  or  a  plot  in  the  form  of  a  bar  diagram 
with  either  a  linear  or  logarithmic  base.  A  logarithmic  scale  is  preferable 
if  the  sound  under  analysis  contains  widely  spaced  components.  If  it  contains 
a  number  of  components  within  a  relatively  selective  (narrow)  frequency  range, 
a  linear  scale  is  preferable. 
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A  sound  whose  analysis  produces  many  discrete  spectral  lines  or  bars 
is  often  called  a  "noise  with  a  line  spectrum"  as  opposed  to  a  pure  tone  of 
single  frequency  and  one  spectral  line. 

4.1.4  Complex  Noise 

Sounds  with  continuous  spectra  (that  is,  noises  which  display  no 
prominent  components  at  discrete  frequencies)  are  generally  analyzed  with 
either  constant  bandwidth  or  constant  percentage  bandwidth  (constant  Q) 
filter  analyzers.  To  be  able  to  compare  the  results  as  obtained  by  different 
analyzers,  the  observed  data  should  be  reduced  to  spectrum  levels.  The 
spectrum  level  is  defined  as  the  level  in  db  of  the  effective  sound  pressure 

in  bands  one  cycle  wide.  This  is  entirely  analogous  to  the  concept  g® /cps 
where  the  mean  squared  acceleration  was  divided  by  the  filter  bandwidth,  to 
reduce  instrumentation  effects. 

To  reduce  a  measured  level  L  db  taken  by  a  filter  band  with  a  width 
of  Af  cps  to  a  spectrum  level  of  S  db,  the  following  operation  is  performed. 


S  =  L  -  10  logjo  Af  db  (128) 


This  spectrum  level  S  is  then  plotted  at  the  geometric  mean  frequency  of  f^ 

equal  to  V f^  fa  of  band  where  f^  and  fg  are  the  nominal  cutoff  frequencies. 

The  function  10  log^^  Af  will  be  a  constant  for  a  heterodyne  type 

(fixed  bandwidth  filter)  analyzer  and  for  a  constant  percentage  bandwidth 
type,  it  will  increase  in  magnitude  by  three  db  for  each  doubling  of  fre¬ 
quency. 


Generally,  all  measurements  of  continuous  spectra  noise  should  be 
reduced  to  spectrum  level. 

4.1.5  Instruments  for  Analyzing  Noise 

In  previous  sections,  various  references  have  been  made  to  filter 
sets,  analyzers,  integrators  and  plots  of  analyzed  data.  Several  types  of  wave¬ 
form  analyzers  and  their  important  characteristics  will  not  be  enumerated; 
some  normally-used  readout  mechanisms  will  also  be  included.  These  instruments 
are  as  follows: 


a.  Octave  Band  Analyzer  In  an  octave  band  analyzer,  both  the 
relationship  between  adjacent  filter  center  frequencies  and  between  the  upper 
and  lower  3-db  points  of  one  particular  filter  is  2  to  1.  (Two  frequencies 


are  an  nth  octave  apart  if  •—  equals  2  .)  A  standardized  set  of  octave 
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passbands,  commonly  used  in  the  United  States  has  the  following  passbands: 

75  to  150,  150  to  300,  300  to  600,  600  to  1200,  1200  to  2400,  and  2400  to 
4800  cps.  Another  commonly  used  series  of  octave  bands  is  50  to  100,  100  to 

200,  _  3200  to  6400  cps.  Besides  the  2  to  1  relationship  between  center 

frequencies  and  cutoff  frequencies,  the  other  useful  relationship  in  connection 
with  octave  band  analysis  is  that  the  3-db  passband  of  any  arbitrary,  constant 
percentage,  octave  filter  is  0.67  times  the  filter  center  frequency. 

b.  Half-Octave  and  Third-Octave  Analyzers  If  a  more  detailed 
analysis  is  required,  a  constant  percentage  bandwidth  analyzer  with  narrower 
fixed  filters  can  be  used.  These  filter  bands  can  be  a  half  or  a  third  of  an 
octave  in  width  or  even  narrower.  The  one-third  octave  analyzers  in  use  are 
normally  designed  on  the  basis  of  the  "preferred  numbers".  A  design  on  the 
preferred  number  basis  selects  center  frequencies  whose  common  logarithm  has 
a  mantissa  which  is  an  even  tenth.  For  example,  .1000,  .2000,  .3000  are  the 
mantissas  respectively  for  1.26,  1.58  and  1.99.  These  antilogs  then  would  be 
rounded  off  to  1.25,  1.6  and  2.0  respectively,  and  they,  as  well  as  the  series 
formed  by  their  multiplication  with  integral  powers  of  ten,  would  form  the 
basis  for  filter  center  frequencies.  Similarly,  those  numbers  whose  common 
log  mantissa  is  .1500,  .2500,  .3500,  etc.,  form  the  basis  for  3-db  crossover 
points  between  filters.  The  third-octave  relationships,  namely  a  factor  of 
1.25  between  adjacent  center  frequencies,  and  a  3-db  bandwidth-to-center- 
frequency  ratio  of  0.231  then  follow  from  the  fact  that  one-third  octave  is 
equal  to  one-tenth  decade. 

Half-octave  relationships  reduce  to  1.414  for  the  adjacent  center 
frequency  ratio  and  0.34  for  the  bandwidth-to-center-frequency  ratio.  (See 
Figure  42.) 

c.  Narrow  Band  Proportional  Bandwidth  Analyzers  Analyzers  which 
have  passbands  even  narrower  than  one-third  octave  are  also  used,  to  determine 
more  exactly  the  presence  of  discrete  components  in  a  signal.  This  type  of 
analyzer  normally  is  continuously  tunable. 

d.  Constant  Bandwidth  Analyzers  Many  different  models  of  constant 
bandwidth  analyzers  are  now  available.  Generally,  all  the  models  are  similar 
in  one  respect,  i.e.,  they  first  combine  the  input  signal  with  a  tunable 
local  oscillator  in  a  modulator  circuit,  and  then  filter  the  result  using  a 
constant  bandwidth  filter.  The  filtering  can  take  place  at  a  fixed  IF  (for 
example,  at  100  kc)  or  at  d-c  (zero  IF).  If  a  fixed  IF  is  employed,  crystal 
filtering  is  often  used,  whereas  in  the  zero  IF  schemes,  Chebychev  or  simple 
resistance-capacitance  (RC)  filtering  normally  is  encountered.  Constant 
bandwidths  as  narrow  as  1/4  cps  operating  to  30  kc  and  higher  have  been 
successfully  used  in  this  type  analyzer.  Other  features  often  offered  in  this 
regard  include:  various  filter  bandwidths,  selectable  by  the  operator; 
optional  detection  of  the  filtered  signal  (linear,  peak,  or  square  law); 
automatic  normalization  (division  by  filter  bandwidths) ;  linear  or  logarithmic 
sweep  of  frequency  versus  time;  mechanical  or  all  electronic  sweep  drives; 
optional  smoothing  of  output  data;  digital  output  of  analyzed  data  and  other 
optional  display  features. 
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on  0=  BANDWIDTH  FACTOR  IN  PERCENT  FREQUENCY 


e.  Parallel  Filter  Analyzer  -  Parallel  Readout  The  constant 
percentage  bandwidth  analyzers  discussed  so  far,  are  examples  of  parallel 
filter  analyzers  which  operate  over  only  one  selectable  passband  at  a  time. 
Other  parallel  filter  analyzers  which  provide  "simultaneous"  readout  from 
many  filters  are  also  available.  One  such  device  uses  480  parallel  magneto- 
strictive  filters  of  constant  bandwidth  and  a  rapidly  commutated  output  which 
can  be  observed  on  an  oscilloscope  or  recorded  on  a  strip  chart. 

Another  method  uses  parallel  constant  Q  filters,  variable  detection 
digital  averaging,  and  a  commutated  output.  These  analyzers  have  the  advantage 
of  presenting  an  almost  continuous  output  over  the  entire  spectrum  of  interest 
and  so  give  an  indication  of  "real  time"  variation  of  the  input  data  spectrum. 

f.  Circulating  Memory.  Analyzers  A  significant  contribution  to  the 
development  of  spectrum  analysis  has  been  achieved  recently  by  the  development 
of  the  circulating  delay  line  analyzer.  At  least  two  different  analyzers 

of  this  general  type  are  presently  available.  One  of  these  analyzers, 
heterodynes  the  incoming  signal  to  a  suitable  frequency  range  where  it  is 
circulated  up  to  360  times  a  second  in  a  crystal  delay  line.  By  a  unique 
phase-advance  arrangement,  a  highly  selective  analysis  is  performed  with  near 
optimum  signal-to-noise  enhancement  and  a  visual  display  of  spectrum  co¬ 
efficients  is  presented  every  second.  In  this  way  spectral  components  over  a 
fixed  range  of  frequencies  are  displayed  with  an  accuracy  of  nearly  one  cps 
in  essentially  real  time. 

Another  analyzer  of  this  type  samples  and  stores  the  incoming  data 
in  several  parallel  delay  lines.  By  using  delay  line,  time -compress ion 
techniques,  a  speed  of  up  to  50,000:  1  is  achieved  and  a  frequency  resolution 
of  0.01  cps,  or  more,  is  available.  A  digital  output  as  well  as  standard 
analog  displays  is  readily  available  from  this  analyzer. 
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Besides  having  the  advantages  of  resolution  and  speed,  these 
analyzers  do  not  present  the  filter  storage  or  ringing  problems  normally 
encountered.  Unless  rapid  observation  of  large  amounts  of  data  is  necessary, 
the  relative  costs  of  this  type  analyzer  might  make  them  impractical. 

4.1.6  Commonly  Used  Readout  Devices 

The  most  commonly  used  readout  devices  in  data  analyzers  are  as 

follows: 

a.  Level  Recorders  High  Speed  level  recorders  are  often  used  to 
display  filter  outputs.  Linear  or  logarithmic  ordinate  scales  are  normally 
available.  With  this  type  of  presentation,  the  time  dependence  of  data 
through  a  fixed  filter  bandwidth  can  be  observed. 

b.  Oscillograph  Recorders  An  oscillograph  record  is  the  most 
straightforward  means  to  display  the  total  unfiltered  signal.  Frequently,  on 
one  record,  are  displayed  not  only  the  total  unfiltered  signals,  but  the 
outputs  of  several  filters  simultaneously  or  the  differentiated  or  integrated 
data  as  well.  Averaging  or  smoothing  or  even  rectification  of  the  filter  out¬ 
put  can  be  done  before  recording. 

c.  Oscilloscope  Oscilloscope  displays  have  been  used  extensively 
with  "real  time"  analysis  systems.  This  has  the  advantage  of  preserving, 
(through  the  use  of  photographic  equipment)  both  the  frequency  and  the  time 
dependence  of  the  data,  since  a  display  of  amplitude  versus  frequency  is  given 
every  display  period.  (Display  periods  presently  range  from  approximately 
every  1/2  second  to  every  two  seconds.  In  some  cases,  continuous  display  of 
frequency -amplitude  information  is  given.)  The  rapid  sampling  and  buffer 
storage  necessary  to  convert  these  visual  displays  to  a  permanent  record 
(other  than  a  photograph)  makes  the  conversion  involved  and  sometimes  expensive 

d.  X-Y  Recorders  In  most  cases  where  a  time  average  presentation 
of  frequency-amplitude  information  is  required,  an  X-Y  plot  is  suitable.  In 
probability  density  analysis,  this  method  has  been  used  almost  exclusively. 

It  can  be  extended  for  use  with  an  oscilloscope  presentation  if  the  process  of 
sampling,  conversion,  and  tabulation  is  previously  undertaken. 

5.0  ANALYSIS  OF  STRAIN  DATA 

5.1  STRAIN  IN  ONE  DIRECTION 

It  recently  has  been  found  that  the  analysis  of  stress  and  strain 
applied  to  structurally  loaded  systems,  can  often  detect  and  sometimes  lead 
directly  to  solution  of  design  or  application  problems  relating  to  structural 
failure.  If  a  simple  structure  is  involved  and  the  direction  in  which  a 
strain  measurement  is  required  is  known,  then  one  gauge  will  give  a  close 
indication  of  true  strain,  "at  a  point".  (The  manner  in  which  strain  is 
measured  with  a  gauge  gives  the  average  over  a  small  distance;  this  normally 
presents  no  problem  unless  the  strain  gradient  is  extremely  large.  If,  then. 
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a  discrete  value  can  be  assigned  to  the  material  under  study  as  its  modulus 
of  elasticity,  stress  can  be  calculated.  Quite  often  however,  because  of 
the  complexity  of  the  structure  to  be  analyzed,  little  or  no  prior  knowledge 
of  directions  of  principle  stresses  is  available.  In  such  cases,  a  strain 
gauge  rosette  is  the  best  approach. 

5.2  STRAIN  GAUGE  ROSETTES 

To  completely  specify  the  state  of  stress  at  any  point  on  a  free 
surface,  three  independent  quantities  must  be  known.  These  quantities  are 
the  magnitude  of  the  two  principle  stresses,  aj  qa  .  and  their  directions 

with  respect  to  some  reference.  Generally,  it  will  be  necessary  to  make  three 
independent  measurements  of  strain  to  determine  these  three  quantities. 

The  three  most  commonly  encountered  strain  rosette  patterns  are 
shown  in  Figure  43.  By  means  of  the  three  or  four  measurements  indicated, 
the  strain  in  any  arbitrary  direction,  with  respect  to  a  fixed  reference,  may 
be  obtained  from  application  of  Mohr's  strain  circle  or  Mohr's  circle  of 
stress . 


The  equations  which  describe  the  principle  stresses,  the  maximum 
shearing  stress,  and  the  direction  of  principle  stresses  in  terms  of  measure¬ 
ment  strains  are  as  follows: 

a.  For  the  rectangular  strain  rosette: 
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Figure  43.  Strain  Rosette  Patterns. 
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where  a,  and  oh  are  the  principle  stress  magnitudes,  t  is  the  maximum 

shearing  stress,  0  is  the  direction  of  the  principle  stresses,  ju  is  Poisson's 

ratio,  E  is  the  modulus  of  elasticity  in  shear  and  e  ,  e,^  ,  e  ,  and  e,  are  the 

abed 

measured  strains. 

Multiple  gauge  connections  can  be  used  to  measure  such  parameters 
as  thrust  on  a  circular  shift,  transverse  shear  on  a  member  in  bending,  torque 
on  a  circular  shaft,  and  bending  moment  or  torque  at  inconvenient  or  in¬ 
accessible  locations. 

5.3  PRESENTATION  OF  STRAIN  DATA 

Dynamic  strain  data,  associated  with  the  structural  evaluation  of 
missile  systems,  can  be  presented  in  various  ways.  For  example,  if  a  launch¬ 
ing  rail  were  instrtimented  with  several  strain  gauges,  the  stress  distribution 
versus  time,  along  the  rail,  could  be  presented  in  tabular  or  pictorial  form. 
Or,  if  an  oscillatory  motion  is  encountered,  as  is  often  the  case,  the  number 
of  cycles  during  which  a  preselected  strain  level  is  encountered  could  be 
presented.  Such  information  would  be  useful  in  fatigue  analysis. 

Since  a  spectrum  analysis  of  strain  data  normally  is  not  required, 
oscillograph  recordings  of  strain  information  in  most  cases  are  sufficient. 
From  these  recordings,  such  values  as  maximum  strain  and  short  term  average 
strain  can  be  determined  as  well  as  correlation  with  other  significant  events. 
The  recorded  data  can  be  sampled  and  converted  into  discrete  numbers.  In  this 
manner,  general  purpose  digital  computers  have  been  successfully  programmed 
to  evaluate  such  parameters  as  those  in  equations  (139)  to  (140). 
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1.  Acceleration  of  Gravity  (g) :  The  acceleration  of  gravity  is  taken  as 

32.2  ft/sec®,  386  in. /sec®,  or  980.7  cm/sec®.  The  nondimens ional  unit  of 
acceleration,  G,  is  defined  as  a/g,  which  equals  the  nimber  of  gravity 
units.  The  acceleration  is  a. 

2.  Amplitude:  In  vibration  terminology,  the  word  amplitude  is  usually 
reserved  to  represent  the  maximirai  value  of  a  sinusoidal  vibration.  If 
used  in  other  senses,  its  meaning  should  be  carefully  indicated. 

3.  Bandwidth,  or  Passband  Width:  Electrical  filters  can  be  constructed  to 
pass  signals  between  two  frequency  limits  and  to  reject  signals  outside 
these  limits.  The  difference  between  these  frequencies  is  the  bandwidth, 
the  average  of  these  two  frequencies  is  the  center  frequency.  Octave 
bands,  half-octave  bands,  etc.,  are  commonly  used.  Smaller  bandwidths 
are  required  if  good  spectral  resolution  is  wanted. 

4.  Conversion  from  One  Band  Pressure  Level  to  Another:  If  the  sound 
pressure  level  is  determined  as  p^  for  a  bandwidth  Afj  >  then  the  band 

pressure  level  for  a  bandwidth  Afg  can  be  calculated  as: 


Afj 

=  Lp^  -  10  log  ^ 


It  has  been  assumed  for  the  above  that  the  spectral  level  is  constant 
throughout  the  frequency  band . 

5.  Decibel  (db)  :  The  decibel  is  a  unit  of  level  which  determines  the  ratio 
of  two  quantities  that  are  proportional  to  power;  the  number  of  decibels 
is  equal  to  10  times  the  logarithm  to  the  base  10  of  this  ratio. 

6.  Discrete  Frequency:  A  vibration  of  discrete  frequency  is  a  sinusoid. 

7.  Distribution  Curve:  The  ordinate  of  the  distribution  is  called  probability 
density  and  the  abscissa  is  the  deviation  of  the  magnitude  of  the  function 
from  its  mean  value.  If  the  function  represents  a  vibration,  and  if  the 
mean  value  of  the  vibration  is  zero,  the  deviation  is  equal  to  the  in¬ 
stantaneous  amplitude  of  the  vibration.  The  area  under  the  distribution 
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curve,  and  bounded  by  two  abscissa  values,  is  equal  to  the  probability 
that  the  instantaneous  amplitude  of  the  vibration  will  be  contained 
between  the  two  abscissa  values.  This  area  can  also  be  considered  as 
the  fractional  part  of  the  time  that  the  instantaneous  amplitudes  of  a 
vibration  are  within  this  range. 

8.  Double  Amplitude:  Double  amplitude  may  be  used  for  sinusoidal  vibrations 
as  equal  to  two  times  the  amplitude.  This  term  has  come  into  use  because 
it  is  the  excursion  visually  observed  in  low-frequency  vibrations. 

9.  Gaussian  (or  Normal)  Distribution:  The  probability  density  for  the  in¬ 
stantaneous  amplitudes  for  the  usual  type  of  random  vibrations  is 
approximately  given  by  the  Gaussian  (or  normal)  distribution  curve. 

10.  Half-Sine  Shock  Pulse:  A  shock  excitation  for  which  the  acceleration 
time  curve  has  the  shape  of  the  positive  (or  negative)  phase  of  one  cycle 
of  a  sine  wave. 

11.  Intensity  Level:  The  intensity  level  of  a  sound,  in  decibels,  is  10  times 

the  logarithm  to  the  base  10  of  the  ratio  of  the  intensity  I  of  this 

sound  to  the  reference  intensity  I  .  The  reference  intensity  1  shall  be 

o  o 

stated.  A  generally  used  value,  especially  for  air  acoustics  is  10 
watts/cm® . 

Note 

Conventional  sound  pressure  meters  and  sound  level 
meters  do  not  measure  intensity.  Hence,  the  words 
"intensity  level"  should  not  be  applied  to  data 
taken  with  them.  Instead,  the  two  expressions 
"sound  level"  and  "sound  pressure  level"  are  used 
in  this  MTP  as  follows:  (a)  Sound  level  is  applied 
to  data  taken  on  instruments  which  meet  the  speci¬ 
fications  for  sound  level  meters  drawn  up  by  the 
American  Standards  Association.  (b)  Sound  pressure 
level  is  applied  to  data  taken  by  a  sound  pressure 
meter  with  a  flat  response.  In  both  cases  the 
reference  pressure  is  0.0002  dyne/cm® . 

12.  Instantaneous  Amplitude:  The  magnitude  of  a  vibrational  motion  at  any 
given  instant  of  time. 

13.  Level:  In  communication  and  acoustics,  the  level  of  a  quantity  is  the 
logarithm  of  the  ratio  of  that  quantity  to  a  reference  quantity  of  the 
same  kind.  The  base  of  the  logarithm  (usually  10),  the  reference 
quantity,  and  the  kind  of  level  must  be  specified. 

14.  Line  Spectra  (Spectrum):  A  line  spectrum  is  an  assembly  of  discrete 
frequencies  plotted  amplitudes. 
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15.  Loudness  Level:  The  loudness  level  of  a  sound  is  numerically  equal  to 
the  sound  pressure  level  in  decibels  of  the  1000-cycle  pure  tone  which 
is  judged  by  the  listeners  to  be  equivalent  in  loudness.  The  1000-cycle 
comparison  tone  shall  be  considered  as  a  plane  sinusoidal  sound  wave 
coming  from  a  position  directly  in  front  of  the  observer.  The  listening 
is  to  be  done  with  both  ears,  and  the  intensity  level  of  the  1000-cycle 
comparison  tone  is  to  be  measured  in  the  free  progressive  wave.  The 
reference  sound  pressure  is  0.000200  dyne /cm?  .  The  unit  is  the  phon. 

16.  Masking:  The  number  of  decibels  by  which  a  listener's  threshold  of 
audibility  for  a  given  tone  is  raised  by  the  presence  of  another  sound. 

17.  Mean  Square  Acceleration  Spectral  Density:  Also  known  as  power  density, 
acceleration  density,  acceleration  spectral  density  or  spectral  density. 
The  term  "acceleration  density"  is  most  frequently  used.  (Velocity 
displacement,  or  other  appropriate  terms  may  be  substituted  for 
acceleration).  As  defined  for  random  vibration,  the  mean  square  accelera¬ 
tion  spectral  density,  or  the  acceleration  density,  is  the  mean  of  the 
square  of  the  acceleration  amplitude  per  unit  bandpass  width;  i.e.,  it 

is  the  limiting  value  of  the  mean  square  acceleration  passed  by  a 
rectangular  bandpass  filter,  divided  by  the  bandwidth,  as  the  bandwidth 
approaches  zero. 

18.  Mean  Square  Acceleration  Density  Spectrum:  Spectral  density  expressed, 
or  plotted,  as  a  function  of  frequency  constitutes  a  density  spectrum. 

The  mean  square  acceleration  density  spectrum  is  abbreviated  to  power 
spectrum,  density  spectrum,  or  acceleration  density  spectrum.  The  last 
term  appears  to  have  most  general  acceptance. 

19.  Microbar  (dyne  per  square  centimeter):  The  unit  of  sound  pressure;  one 
microbar  equals  1  dyne/citr®  . 

20.  Noise:  In  acoustics,  noise  is  defined  as  unwanted  sound.  In  vibration 
terminology,  noise  refers  to  vibrations  which  are  nonperiodic  and  of  a 
generally  random  nature.  This  latter  concept  also  applies  to  background 
electrical  disturbances  that  are  characterized  as  noise. 

21.  Octave:  The  interval  between  any  two  tones  whose  frequency  ratio  is  2:1. 

22.  Peak -to -Peak  Amplitude:  While  the  term  usually  is  S3nionymous  with  double 
amplitude,  it  may  represent  a  different  value  for  periodic  complex  waves 
where  it  represents  the  sum  of  the  absolute  values  of  the  maximum  positive 
and  negative  magnitudes. 

23.  Phon:  The  unit  for  measuring  the  loudness  level  of  a  tone.  The  number 
of  phons  is  equal  to  the  number  of  decibels  a  1000-cycle  tone  is  above 
the  reference  sound  pressure  when  judged  equal  in  loudness  to  the  tone 
in  question. 

24.  Power  Spectrum  Level:  The  power  spectrum  level  is  equivalent  to  the 
spectral  density  described  for  mechanical  vibrations.  It  is  the  acoustic 
power  per  unit  bandwidth.  This  is  proportional  to  the  pressure -squared 
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L  =  L  -10  log  A  f. 
OJS  LL 


25.  Pressure  Spectrum  Level:  The  effective  sound  pressure  for  the  sound 
power  contained  in  a  unit  bandwidth  is  the  pressure  spectrum  level.  If 
A  fg  equals  unity,  (see  equation  in  definition  for  "Conversion  from  One 

Band  Pressure  Level  to  Another),  the  value  of  this  equation  satisfies 
this  definition  and  pressure  spectrum  level  is  seen  to  be: 

Lp^  =  Lp^  -10  log  A 

26.  Probability:  Probability  is  the  likelihood  of  occurrence  of  a  parti¬ 
cular  event.  It  is  generally  estimated  as  the  ratio  of  the  number  of 
occurrences  of  the  particular  event  to  the  total  number  of  occurrences 
of  all  types  of  events  considered.  A  certain  event  has  probability  one 
and  an  impossible  event  has  zero  probability.  The  reverse  of  these 
statements  is  not  necessarily  true,  i.e.,  an  event  probability  one  need 
not  necessarily  happen.  (See  Reference  34.) 

27.  Probability  Density:  Probability  density  is  the  ratio  of  the  change  of 
probability  to  the  change  in  number  of  particular  events  considered. 

In  random  vibration,  this  is  the  ratio  of  the  probability  that  the 
(instantaneous)  amplitude  will  be  within  a  given  incremental  range  to 
the  size  of  the  increment  as  the  increment  approaches  zero  at  a 
specified  amplitude  position. 

28.  Random  Vibration:  A  random  vibration  is  an  oscillation  whose  instantaneous 
amplitude  cannot  be  specified  for  anyone  given  instant  of  time. 

29.  Root  Mean  Square  (rms) :  The  rms  value  of  a  set  of  numbers  is  equal  to 
the  square  root  of  the  average  of  the  sum  of  the  squares  of  these  numbers. 
The  rms  value  of  a  function  f(t)  between  t^  and  tg  is: 


rms 


dt 


J 


30.  Shock  Spectrum:  A  displacement  shock  spectrum  is  the  distribution  of  the 
maximum  relative  displacement  responses  of  a  series  of  single  degree  of 
freedom  systems  to  a  shock  excitation  as  a  function  of  the  frequencies  of 
the  systems.  Shock  spectra  can  be  expressed  in  terms  of  acceleration  or 
velocity  by  multiplying  the  displacement  spectra  by  (2TTf)*  or  2i~f ,  re¬ 
spectively.  Spectra  can  be  obtained  for  either  damped  or  undamped 
systems.  An  undamped  system  is  assumed  unless  otherwise  specified. 
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31.  Sone:  A  unit  of  loudness.  It  is  defined  as  the  loudness  of  a  1000- 
cycle  tone  40  db  above  threshold.  A  millisone  is  one -thousandth  of  a 
sone  and  is  often  called  the  loudness  unit. 

32.  Sound  Intensity:  The  flow  of  a  sound  wave  is  a  passage  of  sound  energy. 
The  power  transmitted  per  unit  of  area  in  the  direction  of  travel  is 
defined  as  the  sound  intensity  and  is  equal  to: 


in  a  free  homogenous  medium  where  I  is  the  density,  C  is  the  velocity 
of  sound,  and  p  is  the  pressure. 

33.  Sound  Power:  If  a  sound  source  emits  an  amount  of  acoustic  power  W, 
then  for  a  free  field  the  sound  intensity  at  any  point  a  distance  r 
from  the  source  is : 

W 

Airr* 

34.  Sound  Power  Level:  The  sound  power  level  of  a  sound  source,  in  db,  is 
10  times  the  logarithm  to  the  base  10  of  the  ratio  for  the  sound  power 
radiated  by  the  source  to  a  reference  power.  The  reference  power  is 
usually  taken  as  one  micromicrowatt,  or  10  watt. 

35.  Sound  Pressure  Level:  The  sound  pressure  level  in  db ,  of  a  sound  is  20 

times  the  logarithm  to  the  base  10  of  the  ratio  of  the  pressure  p  of 

this  sound  to  the  reference  pressure  p  .  The  value  of  p  .  should 

ref  ref 

always  be  stated.  A  common  reference  pressure  used  in  connection  with 
hearing  and  the  specification  of  noise  is  0.000200  dyne/cm^  (or  0.0002 
microbar).  Another  commonly  used  reference  pressure  is  1  dyne/cm®. 

The  two  differ  exactly  74  db .  It  is  to  be  noted  that  in  many  sound 
fields  the  sound  pressure  ratios  are  not  proportional  to  the  square 
root  of  corresponding  power  (intensity)  ratios  and,  hence,  cannot  be 
expressed  in  db  in  the  strict  sense;  however,  it  is  common  practice  to 
extend  the  use  of  the  db  to  these  cases. 

36.  Sound  Spectrum  Level:  The  spectrum  level  is  a  plot,  as  a  function  of 
frequency,  of  20  times  the  logarithm  to  the  base  10  of  the  ratio  of 
the  effective  pressure.  The  unit  is  the  db . 

37.  Spectrum:  A  spectrum  of  a  function  of  time  is  a  resolution  of  its 
components  into  sinusoidal  functions  and  is  expressed  as  a  function  of 
frequency.  The  spectrum  is  also  used  to  signify  a  continuous  frequency 
range  of  sinusoidal  components. 

38.  Standard  Deviation:  The  rms  value  of  the  deviations  of  a  function,  or 
of  a  set  of  numbers,  from  a  mean  value  is  the  standard  deviation.  It  is 
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identified  by  the  Greek  letter  q. 

39.  Sweep  Rate,  Uniform  Sweep  Rate,  Logarithmic  Sweep  Rate,  Constant 
Octave  Sweep  Rate:  A  uniform  sweep  rate  exists  if  the  rate  of  change 
of  frequency  or  df/dt  is  constant.  A  logarithmic  sweep  rate  exists 

if  the  change  of  frequency  per  unit  of  frequency  or  df/fdt  is  constant. 
In  this  latter  case,  the  time  of  sweep  between  any  two  frequencies  of 
fixed  ratio  is  constant.  A  constant  octave  sweep  rate  is  logarithmic 
sweep  rate. 

40.  White  Noise:  White  noise  is  noise  of  random  nature  having  equal  energy 
per  unit  frequency  bandwidth  over  a  specified  total  frequency  band. 
White  random  vibration  is  synonymous  with  white  noise.  A  white  random 
vibration  exists  if  the  mean  square  acceleration  density  curve  is  of 
constant  value  over  the  specified  total  frequency  band. 


MTP  5-1-025 
10  June  1968 


REFERENCES 


1.  Morrow,  C.  T.  and  Muchmore,  R.  B.,  Shortcomings  of  Present  Methods  of 
Measuring  and  Simulating  Vibration  Environments,  Journal  of  Applied 
Mechanics,  Sept.  1955,  pp  367-371. 

2.  Bennet,  W.  R.,  Methods  of  Solving  Noise  Problems,  Proceedings  of  the 
IRE,  May  1956,  pp  609-638. 

3.  Cheng,  D.  K.,  Analysis  of  Linear  Systems,  Addison  Wesley  Publishing  Co., 
1959. 

4.  Schwartz,  M.,  Information  Transmission,  Modulation  and  Noise,  McGraw- 
Hill  Publishing  Co.,  1959 . 

5.  Harris,  C.  M.  and  Crede,  C.  E.  (edited  by)  Shock  and  Vibration  Handbook 
(Volume  II),  McGraw-Hill  Book  Co.,  1961. 


6.  Crandall,  S.  H.,  Random  Vibration,  Technology  Press  of  MIT,  1958. 

7.  Spindler,  C.  W.,  Jr.,  A  Realtime  Complex  Wave  Analyzer,  IRE  -  Proceedings 
of  the  1958  National  Symposium  on  Telemetering,  Section  4.1. 

8.  Morrow,  C.  and  Risen,  D.,  Shock  Spectrvim  Computer  for  Frequencies  up  to 
2000  cps ,  Journal  of  Acoustical  Society  of  America,  January  1956. 

9.  Morrow,  C.  T.,  Averaging  Time  and  Data  Reduction  Time  for  Random  Vibration 
Spectra ,  Journal  of  the  Acoustical  Society  of  America,  Vol.  30,  No‘s  5  & 

6,  June  1958. 

10.  Lee,  Y.  W.,  Statistical  Theory  of  Communication,  John  Wiley  &  Sons, 
Published  1960. 

11.  Harris,  C.  M.  (edited  by)  Handbook  of  Noise  Control,  McGraw-Hill,  1957. 


CM 

Beranek , 

L. 

L. 

,  Acoustic  Measurements, 

John 

Wiley  6c  Sons,  1949 

13. 

Beranek , 

L. 

L. 

,  Acoustics,  McGraw-Hill 

Book 

Co.,  1954. 

14.  Lee,  G.  H. ,  An  Introduction  to  Experimental  Stress  Analysis,  John  Wiley 
&  Sons,  Publisher,  1950. 

15.  Bruel  and  Kjaer,  Audio  Frequency  Spectrometer  Type  2110,  Instructions  and 
Applications  (Manual),  Reprinted  Sept  1959. 

16.  Noise  Measurements  and  Analyses,  Bruel  and  Kjaer  Technical  Review,  No.  1, 
1956. 

17.  Sound  Analysis,  Bruel  and  Kjaer  Technical  Review,  No.  2.,  1957. 


-1- 


MTP  5-1-025 
10  June  1968 


18.  Murray,  W.  and  Stein,  P.,  Strain  Gage  Techniques,  Copywright  by  the 
authors,  1960,  MIT. 

19.  Vigness,  I.,  The  Fundamental  Nature  of  Shock  and  Vibration,  Electrical 
Manufacturing,  June  1959. 

20.  Truxal,  J.  G.,  Control  System  Synthesis,  McGraw-Hill  Book  Co. 

21.  Woodward,  P.  M. ,  Probability  and  Information  Theory  with  Applications 
to  Radar,  McGraw-Hill,  (Pergamon  Press)  1953. 

22.  Grafton,  P.  A.,  Shock  and  Vibration  in  Linear  Systems,  Harper  and  Bros. 
Publisher,  New  York. 

23.  Peirce  and  Foster,  A  Short  Table  of  Integrals,  Fourth  Edition,  Ginn 
and  Co.,  New  York. 

24.  Ratz,  A.  G.,  Ph.  D.,  A  Summary  of  the  Properties  of  the  Rayleight  Distri¬ 
bution  ,  Engineering  Report  #5972,  Ortholog  Division,  Gulton  Industries, 

Inc. 

25.  Summary  Report,  Probability  Density  Analyzer,  Report  No.  895,  Bolt 
Beranek  and  Newman,  Inc.,  Cambridge,  Massachusetts,  May  1962. 

26.  Ratz,  A.  G.,  Ph .  D.,  The  Design  of  an  Amplitude  Probability  Density 
Analyzer  for  Vibration  and  Acoustic  Studies,  Bulletin  0-55,  Ortholog 
Division,  Gulton  Industries,  Trenton,  New  Jersey. 

27.  Landee,  R.  W.,  Davis,  D.  C.,  Albrecht,  A.  P.,  Electronic  Designers  Hand¬ 
book  ,  McGraw-Hill,  1957. 

28.  Stevens,  S.  S.,  Calculation  of  the  Loudness  of  Complex  Noise.  Journal 
of  the  Acoustical  Society  of  America,  Vol.  28,  No.  5.,  Sept.  1956, 

p  807ff. 

29.  Miller,  K.  S.,  Engineering  Mathematics,  Rinehart  Publishing  Co.,  1956. 

30.  Bickel,  H.  J.,  Spectrum  Analysis  with  Delay  Line  Filters,  IRE  Wescon 
Convention  Record,  (Part  8),  1959. 

31.  McIntosh,  V.  C.,  Review  of  Techniques  for  the  Measurement  and  Analysis 

of  Missile  Vibration  Data,  Shock  &  Vibration  Bulletin,  WADC  #24,  Feb.  1957. 

32.  Loudness  Evaluation,  Bruel  and  Kjaer  Technical  Review  #2,  1962. 

33.  Bennett,  W.  R.,  Electrical  Noise,  McGraw-Hill  Book  Co.,  1960. 

34.  Davenport,  W.  B.,  Jr.  and  Root,  W.  L.,  Random  Signals  and  Noise. 
McGraw-Hill  Book  Co.,  1958. 


-2- 


MTP  5-1-025 
10  June  1968 


Bendat,  J.  S.,  Ph .  D.,  Principles  and  Applications  of  Random  Noise  Theory, 
John  Wiley  &  Sons,  Inc.,  1958. 

Graham  D.  and  McRuer ,  D.,  Analysis  of  Non-Linear  Control  Systems, 

John  Wiley  &  Sons,  Inc.,  1961. 

Freeman,  J.  J.,  Principles  of  Noise,  John  Wiley  &  Sons,  Inc.,  1958. 

Morrow,  C.  T.,  Shock  Spectrum  as  a  Criterion  of  Severity  of  Shock 
Impulses ,  Journal  of  the  Acoustical  Society  of  America,  Vol.  29.,  No.  5., 
May  1957. 

Reference  Data  for  Radio  Engineers,  Federal  Telephone  and  Radio  Corp., 

New  York,  N.  Y. 

Kerchner  and  Corcoran,  Alternating  Current  Circuits,  John  Wiley  and  Sons. 

Landee,  Davis  and  Albrecht,  Electronic  Designers  Handbook,  McGraw-Hill 
Book  Company. 

Morrow,  C.  T.  and  Sargent,  H.  1.,  Sawtooth  Shock  as  a  Component  Test, 
Journal  of  the  Acoustical  Society  of  America,  Vol.  28,  No.  5,  Sept.  1956. 

Davison,  S.  and  Adams,  E.  J.,  A  Theoretical  Study  of  the  Multi-Freauencv 
Reed  Gage  for  Measuring  Shock  Motion,  The  David  W.  Taylor  Model  Basin, 
Report  #613,  July  1949. 

Seely,  F.  B.  and  Smith,  J.  0.,  Resistance  of  Materials,  John  Wiley  and 
Sons . 


Strain  Gauge  Measurements,  Bruel  and  Kjaer  Technical  Review  #1, 
January  1957. 


-3- 


